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ABSTRACT 



Performance data from the NASA-Armv OH-6A higher harmonic control 
(HHC) flight test program showed significant reductions in main rotor shaft torque 
and engine power in hover and forward flight [1], The unsteady aerodynamics with 
the higher harmonic control system application, including wake effects, were con- 
sidered to study whether such power reductions are feasible. 

An airfoil oscillating in pure plunge can achieve propulsive force ("Katzinavr 
effect"), as in the case of birds’ wings flapping in flight. Here it will be shown 
that this effect can be enhanced in the presence of layers of shed vorticity with 
the proper phasing. In addition, while it is known that an airfoil oscillating in 
pitch, will typically produce drag at most values of reduced frequency, it is found 
that the presence of another layer of shed vorticity of the proper phase, can 
reduce the drag on the pitching airfoil depending upon wake spacing, reduced fre- 
quency, and phase. Under some conditions the added layer or layers of shed vor- 
ticity will even result in propulsive force acting on the pitching airfoil similar to 
the "Katzmayr" effect for the plunging case. 

It was found, for the OH-6A helicopter, that the measured reductions in 
main rotor shaft torque and engine power are feasible when evaluated with respect 
to the "Katzmayr effect" and the additional drag reduction or propulsive force 
obtained due to pitch and plunge oscillations with the effect of adjacent wake 
layers of shed vorticity. 



/ 

TABLE OF CONTENTS 



I. INTRODUCTION 

A. GENERAL 1 

B. ROTARY WING UNSTEADY AERODYNAMICS 2 

C. ADDITION OF WAKE EFFECT 5 

D. HIGHER HARMONIC CONTROL ( HHC ) 6 

E. OH-6A TEST PROGRAM DATA 7 

F. SCOPE 12 

II. BACKGROUND 15 

III. DRAG REDUCTION DUE TO WING FLAPPING 
( KATZMAYR EFFECT ) 

A. INTRODUCTION 31 



IV 



if. the physics or flapping airfoils in 

LOW-SPEED FLIGHT 32 

C. PANEL CODE EOR STEADY INCOMPRESSIBLE 

IN VISCID FLOW 34 

I). UNSTEADY PANEL CODE 42 

E. TIIEODORSEN'S ANALYSIS 46 

E. GARRICK'S PROPULSIVE FORCE ANALYSIS 47 

G. COMPARISON OF RESULTS 48 

II. ESTIMATED POWER REDUCTION DUE TO IIIIC 

EOR OII-6 A HELICOPTER 52 

IV. ANAIASIS OF WAKE INTERFERENCE EFFECTS 

A. INTRODUCTION 60 

If. UNSTEADY TWO-FOIL PANEL CODE 61 

C. LOEWY'S ANALYSIS 63 

I). EVALUATION OF TWO-FOIL PANEL CODE 68 



v 



E. ANALYSIS OF OSCILLATORY BLADE 



INTERFERENCE EFFECTS 74 

V. PROPULSIVE EFFICIENCY AND FLUTTER ANALYSIS 

A. INTRODUCTION 106 

B. PROPULSIVE EFFICIENCY 107 

C. SINGLE-DEGREE-OF-FREEDOM AIRFOIL FLUTTER 115 

VI. RESULTS AND CONCLUSIONS 

A. DISCUSSION OF RESULTS 120 

B. CONCLUSIONS 132 

C. RECOMMENDATIONS 134 

List of References 135 

Initial Distribution List 142 



VI 



LIST OF TABLES 



Page 

Table (3 1) The OH-6A helicopter rotor parameters 53 

Table (3.2) Cd vs k Tor NACA0015 plunging with h n = 1 " 55 

Table (3.3) Cdav versus reduced frequency kp at hbar = 1 " 59 

Table (4.1) Correction factor for a biplane 72 

Table (6.1) Aerodynamic coefficients vs non-dimensional time (from 

time step 113 to 159) for the case of wake phasing, m =75 
degrees, wake spacing h = 2, reduced frequency k =0.0617 
and plunging amplitude hbar = 0.14 ) 129 

vii 



LIST OF ILLUSTRATIONS 



Fittu. . 1.1) 

Figih. 

Figure (1.3) 

Figure (1.4) 

Figure (2.1) 
Figure (2.2) 
Figure (2.3) 
Figure (3.1) 
Figure (3.2) 
Figure (3.3) 
Figure (3.4) 
Figure (3.5) 
Figure (3.6) 
Figure (3.7) 
Figure (3.8) 



L / Lo phase diagram ( Imaginary vs Real ) 

Power reduction results from higher harmonic control 
open loop test program (NAS A/ Armv/Hughes - Sept. 1982) 
Comparison of 4/rev vertical vibration at the pilot seat versus 
airspeed with HHC "on" and "off”. 

Comparison of 4/rev. lateral vibration at the pilot seat versus 
airspeed with HHC "on" and "off’. 

Loewy's schematic elements of unsteady rotor flow field 
Schematic representation of unsteady rotor flow field. 
Two-dimensional models of unsteady rotor flow. 

Propulsive force on plunging airfoil 

Panel Method Geometry 

Relationship between geometrical quantities 

Panel method representation for unsteady flow 

Comparison of panel code results with Garrick ( k = 0. 1 ) 

Comparison of panel code results with Garrick ( k = 1 ) 

Comparison of panel code results with Garrick ( k = 2 ) 

Blade deflection vs radial station at 2P 



viii 



Figure (3.10) 


Blade deflection vs radial station at 4P 


Figure ( 3.11) 


Blade deflection vs radial station at 5P. 


Figure ( 3 12) 


Cdav vs kp ( plunging amplitude 1" ) 


Figure (4.1) 


Aerodynamic unsteady model for a multi-bladed rotor 


Figure (4.2 a) 


C1/C1 0 vs h/c ( Pane Code ) 


Figure (42b) 


Cl/Cl,, vs h/c ( From Ref. [44] ) 


Figure (4.3 a) 


Local velocity distribution. V 0 /V vs X/C ( Panel Code ) 


Figure (4.3 b) 


Local velocity distribution, V 0 /V vs X/C (after Wegley [46]). 


Figure (4 4 ) 


Correction factor. B vs the chord-gap ratio h/c. 


Figure (4 5 a) 


Pressure distribution vs X/C ( Panel Code ) 


Figure (4 5 b) 


Pressure distribution vs X/C ( Panel Code ) 


Figure (4.6 ) 


Wake pattern (airfoil # 2 plunging at high reduced frequency 



and close to airfoil # 1). 

Figure (4.7) Wake interference between two plunging airfoils in closer 
proximity. 



IX 



Figure (4 8) 
Figure (4 9) 
Figure (4 10) 
Figure (4. 1 1 ) 
Figure (4.12) 
Figure (4. 13) 
Figure (4. 14) 
:re(4 15) 
Figure (4. 16) 
Figure (4.17) 
Figure (4.18) 
Figure (4.19) 
Figure (4.20) 
Figure (4.21) 
Figure (4.22) 
Figure (4233) 
Figure (4.24) 
Figure (4.25) 
Figure (4.26) 
Figure (4.27) 
Figure (4.28) 
Figure (4.29) 
Figure (4.30) 
Figure (4.31) 



Cl vs T* ( m= 0, h = 200) 

Cdvs T* ( m=0. h=200) 

Wake pattern ( m=0. h= 2 ). 

Cl vs T* ( m=0, h=2). 

Cd vs T* ( m=0, h=2 ). 

Wake pattern (m=. 20833, h=2). 

Cl vs T* ( m= .20833 , h =2 ). 

Cd vs T* ( m= 0.20833 , h = 2 ). 

Wake pattern < m = 0.25 . h =2 ) 

Cl vs T* ( m = 0.25 , h=2). 

Cd vs T* ( m = 0.25 , h=2). 

Wake pattern (m = 0.5 , h=2). 

Cl vs T* (m = 0.5 , h=2). 

Cd vs T* (m = 0.5, h=2). 

Average propulsive force coefficient versus wake phasing. 
Wake pattern ( pitch. m = 0, h=2). 

Cl vs T* ( pitch, m = 0 , h =2 ). 

Cd vs T* ( pitch, m = 0 , h =2 ) . 

Cl vs T* ( pitch, m = 0.20833, h=2). 

Cd vs T* ( pitch, m = 0.20833, h=2). 

Cl vs T* ( pitch, m = 0.25, h =2 ). 

Cd vs T* ( pitch, m = 0.25 , h = 2). 

Wake pattern ( pitch. m = 0.5, h =2 ). 

Cl vs T* ( pitch, m = 0.5 , h -2 ) . 



x 



Figure (4.32) 
Figure (4.33) 
Figure (4 34) 
Figure (5.1) 
Figure (5 2) 

Figure (5 3) 

Figure (5.4) 

Figure (5.5) 

Figure (5.6) 

Figure (5.7) 

Figure (5.8) 
Figure (6.1) 



Cd vs T* ( pitch, m =0.5, h =2 ) 

Average drag vs wake phasing in case of pitch 
Average drag vs reduced frequency k. 

Propulsive efficiency q versus 1/k for a plunging airfoil. 
Propulsive efficiency q versus l/'k for an airfoil pitching 
about the leading edge. 

Propulsive efficiency q versus 1/k for an airfoil pitching 
about the quarter-chord point andmid-chord point. 

Propulsive efficiency q versus 1/k for an airfoil 
pitching about the three - quarter- chord point. 

An airfoil restrained to rotate about its leading edge in two- 
dimensional flow. 

Variation with reduced frequency k of the real and imaginary parts 
of the non-dimensional aerodynamic moment my due to pitching 
of an airfoil about its leading edge in an incompressible flow. 

Comparison of Cm (Imag.) vs. k p ( NACA0007 oscillating in 
pitch about the leading edge for ct 0 = 1° and N=200 panels ). 

Imag. (my) vs k (for an airfoil pitching at ea = -0.5. 0.0,0. 5, and 10) 

Drag coefficient vs angle of attack. 



XI 



NOMENCLATURE 



angle of attack 
pitching amplitude 
air density 
rotor angular speed 
frequency of oscillation 
natural frequency of oscillation in plunge 
natural frequency of oscillation in pitch 
coordinate of axis of rotation 
airfoil half-chord 
airfoil chord 
airfoil lift coefficient 
reference airfoil lift coefficient 
second airfoil lift coefficient 



airfoil drag coefficient 



Cd av 

011 

012 
C m 

Cpxh 

Cpa 

D 

h 

h 0 

k 

kp 

L 

L a > L h 

M 

Mh , M a 



average drag coefficient 
reference airfoil drag coefficient 
second airfoil drag coefficient 
moment coefficient 

propulsive force coefficient due to plunge 

propulsive force coefficient due to pitch 

average drag force 

non-dimensional wake spacing 

non-dimensional plunging amplitude 

reduced frequency ( based on airfoil half-chord ) 

reduced frequency ( used by panel code and based on airfoil chord ) 

lift force 

aerodynamic coefficients used for Theodorsen analysis 
moment 

aerodynamic coefficients used for Theodorsen analysis 



xiii 



m ratio of frequency of oscillation to rotor speed 

N number of blades 

Pxh propulsive force due to plunge 

P power 

R rotor radius 

r local radius 

S blade area 

t* non-dimensional time 

U, v free stream velocity 

x 0 coordinate of axis of rotation 



XIV 



Acknowledgments 



This research was conducted at the Department of Aeronautics and Astronautics. 
Naval Postgraduate School, Monterey, California, under supervision of Prof E. 
R. Wood . Chairman of my doctoral committee. Prof. M. F. Platzer, Prof. G. V. 
Hobson. Prof. D. Danielson (Department of Mathematics), and Prof. J. V Sanders 
(Department of Physics). 

I would like to express my sincere appreciation to Prof. E. Roberts Wood 
who provided me with this great opportunity to work in this area and for his 
continuous support, ideas, guidance, knowledge and encouragement. Special thanks 
and deep gratitude are also due to Prof. Max F. Platzer for his continuous support, 
guidance, knowledge and valuable hours of council. Also I would like to express 
my gratitude and appreciation to Prof. Garth Hobson for his help and support, 
guidance, encouragement, and improving my knowledge in classes and in conducting 
this research. 

Deep appreciation goes to Prof. Donald Danielson for his sincere guidance 
and the great pleasure to discuss and understand the physics behind the mathematical 
problems during the hours of his special and useful Applied math, class. Special 
thanks to Prof. James Sanders and Prof. Allan Coppens for their great help, guidance 
and valuable knowledge they provided me in the Department of Physics. 

Special appreciation is due my mother, my brother and my sister. Without moral 
support of my wife, Neven, my lovely daughter, Yasmin, and my nice son Aboda, 
this work could not have been carried out. To them I dedicate this thesis. 



xv 



Summary 



Open loop performance data from the NASA-Army OH-6A higher harmonic 
control (HHC) flight test program. Wood et al.[l], show significant reductions 
in main rotor shaft torque and engine power in the airspeed regime from hover 
to 100 knots. Depending upon HHC controller phase and helicopter airspeed, 
reductions in power were recorded as large as 20%, with reductions of the order 
of 10% being typical. 

Higher harmonic control is an active helicopter vibration control concept 
which alters the aerodynamic loads on the rotor blades such that the blade re- 
sponse is reduced. This in turn reduces the vibratory forces and moments acting 
at the hub, which cause airframe vibration. Basically, HHC is an electronic, 
computer-controlled active vibration suppression system which senses and cancels 
vibrations in a helicopter by N/revolution feathering or pitch motion of the rotor 
blades, N being the number of blades 1 . Very substantial vibration reduction was 
achieved under a NASA-Army sponsored program using a modified OH-6A he- 
licopter. Higher harmonic blade pitch control was implemented by superimposing 
4/rev. (32 Hz) swashplate motion upon basic cyclic and collective control inputs. 



1 E. Roberts Wood, R. W Powers. J. H. Cline, and C E. Hammond. "On Developing and 
Flight Testing a Higher Harmonic Control System." Journal of A.H.S January. 1985. 
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During the program, in addition to reducing \ibration levels, the HHC system 
showed the potential for decreased helicopter power requirements. Recent studv 
of the OH-6A flight test data indicates that the mechanism by which the power 
reduction was achieved is related to the unsteady aerodynamics associated with 
HHC. This is supported by the fact that higher harmonic control by the very 
nature of the method by which it is implemented through oscillating an airfoil 
in pitch with resulting plunge motion, requires definition of unsteady aerodynam- 
ics to adequately model its physics 1 

The fundamental closed form solutions of Theodorsen [3], and Loewv [4] 
provide the basis for theoretical work in this area. The closed form theory 
shows rapid changes in the lift deficiency function with changes in reduced fre- 
quency, wake spacing and frequency ratio. In the past, emphasis in the study 
of unsteady aerodynamics has tended to focus on flutter instability and the effect 
of unsteady aerodynamics on generating lift and torsional loads. 

In this research, we are interested in performance, and emphasis will be on 
the effect of unsteady aerodynamics on the drag of the airfoil. When the drag 
force is reversed, acting to propel the airfoil forward, it is generally referred to 
as either negative drag or propulsive force. The classic reference on this subject 
is that by Garrick [5] 



1 E. R. Wood. Max F. Platzer. Ahmed Abourahma. Mark A. Couch " On the Unsteady 
Aerodynamics of Higher Harmonic Control." Paper No. C 17. Nineteenth European 
Rotorcraft Forum. Cemobbio (Como), Italy, September 1 993. 
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While Garrick’s [5] work shows that an airfoil oscillating in pitch will 
typically produce drag in the lower reduced frequency range, k. it is found that 
the presence of another layer of shed vorticitv of the proper phase can reduce 
the drag on pitching airfoils and depending upon wake spacing, reduced frequen- 
cy, and phase, may even enhance the propulsive force acting on the airfoil similar 
to the "Katzmayr " effect for plunging airfoils. 

The specific mechanisms, that have been used in this research to discuss 
die phenomenon of OH-6A helicopter rotor power reducdon. are: ( 1) the unsteady 
-.odynamics associated with a plunging airfoil and how that produces a propul- 
sive force or decrease in rotor power ; (2) the influence of a layer of wake 
vorticity on this propulsive force, which represents the influence of the unsteady 
aerodynamics of one rotor blade on adjacent rotor blades, (3) the unsteady aero- 
dynamics associated with an airfoil oscillating in pure pitch and how this pro- 
duces in general a drag force on the airfoil unless accompanied by layers of shed 
wake vorticity ; and (4) the influence of the phasing of a layer of shed vorticity 
on this drag force, and how, depending upon specific phasing, the effect on the 
airfoil will be either drag or propulsive force. 

The analysis of the above effects was based upon; (1) the classical flat 
plate unsteady aerodynamic theory for an incompressible fluid of Theodorsen [3], 
its extension to include chordwise forces (drag) by Garrick [5], the rotor wake 
analysis of Loewy [4], (2) the panel method developed bv Platzer et al. [6], 
that can be applied fo r unsteady incompressible flow past airfoils or airfoil com- 
binations of arbitrary geometry and amplitudes of different kinds of motions. 
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In the first approach. Theodorsen's well known lift deficiency function C(k) 
has been modified by Loewv for the case of a hovering rotor. Airfoil drag or 
propulsive force was determined by the methods of Garrick [2], In the second 
approach, the panel code allowed a systematic study of the lift and drag produced 
by oscillating airfoils and airfoil combinations in an incompressible inviscid 
flow. 

The panel code results were checked against different references. The results 
showed good agreement with these references. For the parameters of the OH- 
6A rotor in hover, the main conclusion was that the reductions in drag and 
propulsive forces generated on an airfoil, oscillating in pitch or plunge, due to 
wake phasing, are sufficient to explain the reductions in power recorded on the 
OH-6A during the NASA- Army flight test program [1], 
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I . INTRODUCTION 



A. GENERAL 

The determination of the unsteady aerodynamics associated with a helicopter rotor has 
presented a challenge to the helicopter analyst since the conception of the helicopter. The 
effect of blade oscillations on resulting airloads has prompted new interest in oscillatory 
blade aerodynamics in connection with possible performance gains. Rotary wing unsteady 
aerodynamics is considerably more complex than fixed wing aerodynamics. 

The simplest case is that of hover, but even here the problem is complicated by the 
influence of the shed helicopter wake from the rotor. In forward flight, the problem is 
further complicated by the fact that the velocity of a local blade element ranges at the 
blade tip from transonic on the advancing side to low subsonic on the retreating side, and 
at the blade root from low subsonic on the advancing side to reversed flow on the retreat- 
ing side. 

For this reason there is a great need for advances in the present theory' in order to 
explain modem helicopter problems, yet, the number of advances in rotorcraft unsteady 
aerodynamic theory are very limited, so we need to review current unsteady aerodynamic 
theory for rotorcraft and to advance the theory in one of the areas where improvement is 
badly needed. 

The goal of this thesis is to explain the measured power reductions in the OH-6A 
higher harmonic control (HHC) flight test program conducted by E. R. Wood et al [l]. 
This research will consider the unsteady aerodynamics with HHC system application, 

including wake effects, in order to determine whether such power reductions are feasible. 

1 



B. ROTARY WING UNSTEADY AERODYNAMICS 



In this section, some elements of unsteady airfoil aerodynamics will be discussed. 
Consider an airfoii undergoing sinusoidal pitch and plunge motion. As the airfoil performs 
an oscillation, vortices are shed into the medium with a circulation strength equal in mag- 
nitude to the increase in circulation about the airfoil, but opposite in direction. These dis- 
turbances are stored in the fluid because the shed vorticitv convects downstream at the 
local flow field velocity. 

The counter-rotating vortices induce a sinusoidal How field which further changes the 
.t. The result of this complex flow field is a time difference or delay between the 
airfoil's motion and the induced aerodynamic forces. This delay is known as the phase 
lag, y [34] . To simplify calculations for this type of motion, it is common to describe 
the airfoil position and the associated aerodynamic forces by complex variables. For pure 
plunge oscillations the vertical motion of the airfoil is described by the real part of the 
following equation: 

h(t) = h 0 e ioJt (1.1) 

where h 0 is a complex number, and to is the frequency of oscillation. Similarly, the lift is 
de < "‘ j by : 



L= L 0 r e 



1 cot 



where L 0 is the quasi-steadv lift given by the expression : 



(1.2) 



L 0 = .5 pU 2 S C, a (h / U) 



(1.3) 



2 



This is termed quasi-steady because the angle of attack is represented by h/U. 
The values of r 0 and \\i represent the magnitude and the phase respectively, of the 
true instantaneous lift relative to the quasi-steadv lift. The variables r 0 and \\) in 
general depend on the reduced frequency k. the Mach number M. and the Reynolds 
number. For inviscid. incompressible fluid the values of r 0 and \j will be only a 
function of k [34], A complete solution for the oscillating flat-plate airfoil in an 
incompressible inviscid flow has been obtained by Kussner and Theodorsen and 
reproduced from Reference [35] in Figure (1.1 ) 
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Figure (l.l) 



L/L () Phase Diagram ( Imaginary 



vs Real ) . 
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C. ADDITION OF WAKE EFFECT 

For the helicopter in hover, the rotation of the wing introduces a number of features 
that require special attention, notably the returning vortex wake, the time varying free 
stream, and radial flow, and a fundamentally transcendental geometry that requires either 
approximate or numerical solutions. The lift on a wing is due to its bound circulation. 
Conservation of vorticitv in the tlow requires that there be trailed and shed vorticity in 
the wake behind the wing. 

The spanwise variation of bound circulation results in trailed vorticity parallel to the 
free stream direction. Time variadon of the bound circulation leads to shed vorticity par- 
allel to the wing span The wake is composed of sheets of vorticity convected downstream 
from the trailing edge by the free stream velocity. The wake of a rotor in hover or in ver- 
tical flight consists of helical vortex sheets below the disk, one from each blade. Unsteady 
motion of the rotor blade will produce a shed vorticity in the wake spirals. With low disk 
loading ( at low collective pitch settings ), the wake remains near the rotor disk and there- 
fore passes close to the following blades. Thus the wake vorticity is not convected down- 
stream of the airfoil as with a fixed wing, and the shed vorticity sheets below the rotor 
disk must be accounted for to correctly estimate the unsteady loads. 

For high inflow (high collective pitch settings ) or in high-speed forward flight, the 
wake will be convected away from the blades. The returning shed wake influence is pri- 
marily a concern of hover and vertical flight, and forward flight at airspeeds less than 
100 knots. Assuming high aspect ratio of the blade, lifting line theory requires a knowl- 
edge of the loads on the blade section, and the returning shed wake of the rotor must be 
incorporated into the two-dimensional unsteady airfoil theory. 
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The wake far from the blade section will have little influence, so the emphasis is on 
modeling the wake near the blade, which for low inflow consists of vortex sheets that are 
nearly planar surfaces parallel to the disk plane. Based on these considerations, a two- 
dimensional model for the unsteady aerodynamics of the rotor can be constructed. Loewv 
in 1957, developed a two-dimensional model for the unsteady aerodynamics of the blade 
in a hovering rotor., including the effect of the helical shed wake. 

D. HIGHER HARMONIC CONTROL (HHC) 

The primary source of helicopter vibration is the higher harmonic blade loads gener- 
ated by the aerodynamic environment at the rotor disk. For an N-bladed rotor, the oscil- 
latory (N-l)P, NP, and (N+1)P blade loads are transmitted as exciting hub forces and 
moments to the airframe at a dominant frequency of N per rev (NP). A significant reduc- 
tion in vibration level has been achieved by industry in the last twenty five years. The 
desired goal was to achieve 0.02 g and this goal cannot be met without a quantum 
advance in vibration control technology. The vibration levels first specified by the U. S 
.Army for the AAH/UTTAS procurement in the early 1970s were later revised upward in 
the mid 1970s to better reflect realistic design goals consistent with the state-of-the-art in 
helicopter vibration control [1]. 

A higher harmonic control system drives the blades in pitch at the (N-l)P, NP, and 
(N+1)P harmonics of rotor rotational speed, generating new unsteady airloads, which in 
combination with the new oscillatory inertial loads, cancel the harmonics of blade loads 
that cause airframe vibration. Therefore the vibrations are suppressed at the source. 
Although there are several ways to implement HHC on a rotor system, the approach 
generally followed to date has been by blade root feathering using swashplate oscillations . 
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By means of electro-hydraulic servo-actuators, the svvashplate is excited in the 
collective, longitudinal cyclic, and lateral cyclic modes at NP resulting in blade pitch 
oscillations at three distinct frequencies of (N-l)P, NP, and (N+1)P in the rotating frame. 
In the system considered, developed for the four-bladed OH-6A helicopter, higher har- 
monic blade feathering for vibration reduction is achieved by superimposing 4/rev swash- 
plate motion upon basic collective and cyclic flight control inputs. 

Perturbing the stationary swashplate at 4/rev both vertically and in pitch and roll 
results in third, fourth, and fifth harmonic blade feathering in the rotating system. Fourth 
harmonic blade feathering is achieved by oscillating the swashplate vertically about its 
collective position, while third and Fifth harmonic blade feathering in the rotating system 
results from 4/rev tilting of the stationary swashplate in pitch and roll about its cyclic tilt 
position [2], Results of experimental efforts by Hammond [52] showed that successful 
suppression of vibration can be achieved by oscillating the blades at relatively small 
angles. 

E. OH-6A (NASA/ARMY/HUGHES) TEST PROGRAM DATA 

The goal of this research has been a systematic study of the OH-6A helicopter 
HHC flight test program, in order to understand and explain the significant reduction in 
main rotor shaft torque and engine power that was measured in flight. Shown in Figure 
( 1.2 ) are the power reduction results from the HHC Open Loop Flight Test Program 
(NASA/ Army/ Hughes - Sept. 1982). 
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Mgiire ( 1.2 ): Power Reduction Results from Higher Harmonic Control (111 1CT) 
Open Loop lest Program (NASA/Anny/lluglies - Sept. 1982) 
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Figuics (1.3) and (1.4) (From Ref f 1 ]) show the veiy substantial vibration reduction 
achieved under a NASA-Armv sponsored program using a modified OH-6A helicopter. Higher 
harmonic blade pitch control was achieved by superimposing 4/rcv. (32 Hz) swashplate 
motion upon basic cyclic and collective control inputs The aircraft was flown from zero 
airspeed to 100 knots with the HHC system operated both open loop (manually) and closed 
loop (computer controlled) 




Figure ( 1.3 ): Comparison of 4/rev Vertical Vibration at the Pilot 
Seat Versus Airspeed with HHC "oil" and "off. 
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Figure ( 1.4 ): Comparison of 4/rev. Lateral Vibration at the Pilot Seat 
Versus Airspeed with HFIC "oil" and "off. 
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The significance ot' power reduction means that: (1) the helicopter payload can be 
increased. (2) the maneuverability can be markedly increased. (3) marginal operating condi- 
tions. such as at sea in high sea states can be better met. and (4) the helicopter operating 
costs can be substantially reduced. Increased interest in HHC and its potential benefits for 
the helicopter now require that we try to use the proper tool to understand this phenome- 
non. It would appear that we presently have the computational toois. namely the panel 
method for unsteady inviscid. incompressible flow past airfoils or airfoil combinations of 
arbitrary geometry, developed by Platzer et al [6], and the single wake dosed form anal- 
ysis developed from Loewv’s theory by Wood and Couch [40], 
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F. SCOPE 



Chapter II contains a quick review of the most important research starting with the 
work of Knoller[10] in 1909 and Betz [l I] in 1912. They explained the bird’s ability to 
generate forward thrust by means of wing flapping. In 1922. R. Katzmavr, experimentally 
verified the work of Knoller and Betz. Birnbaum, von Karman. Theodorsen. and Garrick in 
the years 1924 through 1936. were the pioneers of flapping wing aerodvnamics. 

The fundamental closed form solutions by Theodorsen [3], Loewv [4], Garrick [5], pro- 
vide our basis for theoretical work in this field. While Loewv s work on wake-induced 
flutter helps to explain the phenomenon, it also points up difficulties to be overcome. That 
is the closed form theory shows rapid changes in the lift deficiency function with changes 
in reduced frequency, k; wake spacing, h: and frequency ratio, m. 

Garrick further developed the work of Theodorsen by deriving the horizontal force for- 
mulas for airfoils or airfoil-aileron combinations oscillating in any of the three degrees of 
freedom: vertical flapping, torsional oscillation about a fixed axis, and angular oscillation 
of the aileron about a hinge. Also presented in Chapter II. Loewy's two-dimensional model 
for the unsteady aerodynamics of the blade of a hovering rotor, is presented to provide the 
base for discussing the obtained results. Also presented in Chapter II. is a quick survey of 
the most relevant references who developed and used panel methods to study unsteady 
aerodynamics of oscillating airfoils in an incompressible, inviscid flow. 
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In the third chapter, the propulsive force is estimated according to the OH-6A helicop- 
ter parameters obtained from reference [36] In this analysis, the blade is divided into sec- 

tions. Assuming hovering conditions, the reduced frequency at each section was evaluated. 
Using the unsteady panel code the time history of the airfoil drag ( or propulsive force ) 
coefficient was recorded. For the assumed plunging motion with plunging amplitude as 
obtained from reference [1], for a given mode of vibration, the average propulsive force 
for each section was calculated. 

The total average propulsive force was obtained by adding together the contribution 
from each blade section. The obtained results showed that the resulting propulsive force 
from a single airfoil oscillating in plunge is not enough to explain the significant power 
reduction on the OH-6A helicopter. This led us to study the aerodynamic interference 
between two blades, which is the main subject of the fourth chapter. 

In Chapter IV, the two-airfoil unsteady panel code developed by Pang [ 25 ], is used 

to study the wake interference between two airfoils. This code has been formulated to 

solve for the potential flow for two airfoils executing unsteady motion in an inviscid 
incompressible flow medium. It is an extension of the code, developed by Teng [24], for 
single airfoils. The technique uses the well known panel method for steady flow and 
extends it to unsteady flow by accounting for the continuous shedding of vortices into the 
trailing wake. 
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Numerous case-runs are presented to illustrate the aerodynamic interference between two 
oscillating airfoils. Results obtained for plunging and pitching airfoils show a significant 
production of propulsive force at certain wake phase angles. In Chapter V the question of 
propulsive efficiency due to plunging or pitching and the possibility of single-degree-of- 
freedom pitching flutter are discussed. Chapter VI discusses the findings of this research 
and presents the main conclusions. 
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II. BACKGROUND 



The effect of plunge and pitch oscillations on the resultant drag of an airfoil 
has been the object of aeroelastic research back almost to the time of the Wright 
brothers (1903) In 1909 and 1912. respectively, Knoller [10], and Betz [11] were 
the first ones to explain the birds' ability to generate a forward thrust by means 
of wing flapping 1 . Another one of the earliest investigators of this era was R. 
Katzmavr [7], who in 1922 conducted two series of experiments to explore the 
effects of oscillatory fluid or airfoil motion on resultant lift and drag. 

Katzmayr's experiments which were performed in the aerodynamic laboratory 
of the Technical University of Vienna, may be divided into two series. In the first 
series, the angle of attack of the wing model was changed by oscillating the wing 
model about an axis parallel to the span and at right angles to the airflow. In the 
second series of experiments, the airfoils were stationary while the airflow itself 
was subjected to periodic oscillations. Here, for certain air velocities and frequen- 
cies, and considering thick airfoil sections. Katzmavr measured negative drag or 
positive propulsive force. This discovery, to be later known as "Katzmavr effect", 
helped to explain the propulsive mechanism of birds flying through the air. 



1 Platzer. M. F.. Neaee. K. S . and Pang, C. K.," Aerodynamic Analysis of Flapping Wing 
Propulsion". AIAA paper No 93-0484. 31st Aerospace Sciences meeting, Reno. NV, January 
11-14, 1993. 
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Later. E. G. Richardson [8] used the"Katzmavr effect" to explain the locomotion 
of fish through water. In his paper, Katzmavr summarized his results to show:(l) 
the effect of flowing air, whose direction is undergoing constant periodical change, 
is quite favorable, (2) the wing section which exhibits favorable characteristics in a 
constant airflow, works better in an oscillating How. 

In 1935, Theodore Theodorsen [3] introduced the general theory of aerody- 
namic instability and the mechanism of flutter. In the first pan of his paper. 
Theodorsen developed the velocity potentials due to flow around oscillating air- 
foil-ailerons. In the later pan of the paper, he developed the differential equations 
of motions and determined the flutter speed. 

The most relevant work with respect to the present research is that by I. E. 
Garrick [5], who applied the work of Theodorsen [3] to obtain a closed form ana- 
lytical expression for the propulsive force generated by an oscillating airfoil. An 
airfoil oscillating in pure plunge or flapping motion was shown by Ganick to have 
a propulsive force throughout the entire range of reduced, k, while for an airfoil 
oscillating in pitch, propulsive force can only be achieved for a specified range of 
reduced frequency, k, where k must be greater than 2.0 . 



16 



In 1942, W Schmidt [14] recognized that the Knoiler-Betz or Katzmavr effect 
applies to both a flapping airfoil in uniform flow or a stationary airfoil in oscillat- 
ing flow. He demonstrated experimentally that a stationary airfoil positioned in the 
wake of a flapping airfoil increases the propulsive efficiency to almost 100 percent 
throughout the whole frequency range because the stationary airfoil converts the 
vortical energy generated by the flapping airfoil into additional thrust. He also 
found that the mechanical limitations inherent in pure flapping motions could be 
overcome by an arrangement which he called the wave propeller. 

H. Bosch [15] investigated the interference of two lifting surfaces in two- 
dimensional unsteady incompressible flow, using classical flat -plate theory. He 
showed that a harmonically flapping airfoil upstream of a stationary airfoil increases 
the propulsive efficiency to almost 100 %, thus confirming Schmidt’s experimental 
findings. Bosch’s analysis is limited to flat-plate airfoils oscillating about a zero 
incidence mean position. 

In 1957, Loewy [4] developed a two-dimensional model for the unsteady 
aerodynamics of the blade of a hovering rotor. He derived an expression for the 
two-dimensional oscillatory lift deficiency function as a function of the reduced 
frequency, wake spacing and the ratio of the airfoil frequency of oscillation and the 
rotational speed. 
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The wake in hover or vertical flight consists of helical vortex sheets below the 
disk, one from each blade. Unsteady motion of the rotor blade will produce shed 
vorticity in the wake spirals. 

In Figures ( 2 . 1 A ) through Figure (2. ID), from Ref. [4], Loewv illustrated the 
basic elements of the unsteady flow field of a helicopter rotor. A rotor in hover- 
ing or vertical flight produces a trailing tip vortex with a downward axial veloc- 
ity that, if otherwise undisturbed, forms a contracting helix as shown in Figure 
(2.1A). If the inflow over the rotor disk is constant, then the fluid off the trailing 
edge of the blades makes a helical surface with horizontal radial elements as shown 
in Figure (2. IB). 

If there is an oscillation in blade effective angle of attack, blade lift will alter- 
nate also, and as a result of these changes in lift, vortices will be shed continu- 
ously at the blade trailing edge. These vortices fall along the horizontal radial 
elements of the helical surface shown in Figure (2. IB), so long as the oscillations 
in angle of attack are small. Figure (2. 1C) illustrates this local sheet of shed vor- 
ticity. It should be noted that vorticity is considered to be on the helical surface 
shown in Figure (2. IB); the vertical displacements from that surface shown in Fig- 
ure (2. 1C) represent the strength of the vorticity at a particular azimuthal and radial 
position [4], 



18 




Figure (2.1) Eoewy's Schematic Elements of Unsteady Rotor Flow Field 
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The variation in this vertical displacement (hence vortex strength) around the 
azimuth corresponds to the history of the motion of a given blade element at a 
fixed radius. Variation of shed vortex strength in the radial direction at any fixed 
azimuth angle is a function of the variation with blade span of (1) blade pianform; 
(2) amplitude of oscillation of effective angle of attack, and (3) relative air veloc- 
ity. Since (from the Helmholtz theorem ) vorticitv cannot begin or end in space, a 
variation of shed vorticitv in the radial direction implies that there are trailing vor- 
tices at constant radii similar to and inboard of the tip vortex. These trailing vorti- 
ces have been induced in Figure (2. ID). 
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The schematic drawings. Figures (2.1 A) through Figure (2. ID), Ref. [4], indicate 
the difficulty of attempting to obtain a complete representation of unsteady rotor 
aerodynamics. One means of simplifying the problem is to consider extreme "pitch" 
values for the helix. For the case of very high inflow velocity, u, in relation to Q, 
and the opposite, where u is very low compared to Q This is shown in Figure 
( 2 . 2 ). 

When the vertical spacing between adjacent helical surfaces of shed vorticity is 
very large, then we expect that all shed vorticity beyond a small fraction of a rev- 
olution would be too far below the blade in question to have a significant effect. 
Under these conditions, it would be sufficient to account for only the attached 
vortex sheet within that fraction of a revolution, as in Figure (2. 2A). 

On the other hand, when the vertical spacing between the adjacent helical 
surfaces of shed vorticity is very small, all the sheets of shed vorticity tend to pile 
up on each other, and the effect of that vorticity close to the blade in question 
( shed by blade passes ) is of more importance than that which exists beyond 
a reference azimuth angle on either side of the blade. This situation is depicted in 
Figure (2.2B), Ref. [4], 
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A. High Inflow 



B. Low Inflow 



figure (2.2) Schematic Representation of Unsteady Rotor Flow Field. 
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With low rotor thaist values ( low collective pitch ) the wake remains near 
the rotor disk and therefore passes close to the following blades. Thus, the shed 
vorticity sheets below the rotor disk must be accounted for to correctly estimate the 
unsteady loads. For high inflow or high-speed forward flight the rotor wake will 
be convected away from the blades. The wake far from the blade section will have 
little influence, so emphasis is on modeling the wake near the blade. In the case 
of low inflow, this consists of vortex sheets that are nearly pianar surfaces paral- 
lel to the disk plane [4] 

Based on these considerations, a two-dimensional model for the unsteady aero- 
dynamics of the rotor can be constructed. It is assumed that the chord, amplitude 
of the oscillation in effective angle of attack, and relative airspeed vary slowly 
enough with span so that what occurs aerodynamically at one blade radius station 
is essentially duplicated on either side of it. 

The assumption means that the flow problem at a given blade section is two- 
dimensional. Furthermore, consistent with the idea that only the vorticity near the 
blade section has an important effect, one may allow the planar rows of vorticity to 
extend to infinity in the horizontal direction in order to achieve mathematical sim- 
plification. It should be stated that this is a first-order theory, so that such effects 
as that of the wake upon itself have been neglected. The oscillatory flow induced 
at the airfoil in the free-steam direction, for example, will affect the airfoil lift and 
moment as well as the vorticity shed at the trailing edge. 
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Since the shed vorticity is being accounted for as perturbations in the main 
How. this component of induced velocity will be small compared to the free-steam 
velocity (except when the spacing approaches zero) and must react with small, first- 
order airfoil displacements. Because of the physical unrealities associated with the 
standard mathematical model representation of the problem (such as zero thickness 
of wake and airfoil, zero displacement of the wake from horizontal planes, and 
infinite velocities at the core of a vortex), the accuracy of results as the wake 
spacing goes to zero is questionable. The effect of the horizontal component of 
induced flow, therefore, has been neglected in Loewv's analysis. 

Shown in Figure (2.3) is the two-dimensionalized model of unsteady rotor flow 
[4], The portion of the circular surface which is determined by (1) a particular 
blade radius, (2) the azimuth angle on either side of the blade section (within 
which the shed vorticity is of importance), and (3) the vertical distance spanned by 
a given numbers of rows of vorticity, can be regarded as a plane, one in which 
the two-dimensional unsteady aerodynamic problem may be attacked. 

This is shown in Figure (2.3 A), with an arrow indicating that the rows of 
vorticity under the rotor disc extend to infinity. Consistent with the idea that only 
the vorticity near the blade section has an important effect, one may allow the pla- 
nar rows of vorticity to extend to infinity in the horizontal direction in order to 
achieve mathematical simplification. These final modifications are shown in Figure 
(2.3 B), which is a two-dimensional model. Ref. [4], of unsteady rotor aerodynamics 
lor a single-bladed rotor operating at low inflow. 
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Figure (2.3) Two-dimensional Models of Unsteady Rotor Flow. 
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The general approach used by Loewv is the same as that of Schwarz [30]. 
Only the effects resulting from oscillatory motion of a thin airfoil are considered, 
and the airfoil itself is represented by continuous vorticity distributed along a 
straight line [9], It may be also noted that the effect of viscosity in damping the 
shed vorticity could be accounted for by considering only a finite number of rows 
beneath the reference airfoil, or by multiplying the vortex strength by an experi- 
mental decay function related to the time through the rotational speed and the 
number of rotational cycles corresponding to a particular vortex row. 
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This has been omitted front the analysis since vortex decay time is long 
compared to the rotational periods of most rotors. An extension to Loewv ’s work is 
that by Shipman and Wood [29] in 1970. The authors present a method for 
determining rotor blade flutter in forward flight. Incompressible unsteady aerody- 
namic theory was applied where shed vorticity both in and below the plane of 
the rotor is accounted for. This was made possible by noting that the damping for 
rotor blade aeroelastic instability varies with velocity Since the tangential velocity 
of any blade segment varies with azimuth, so will flutter damping'. 

It is assumed at the onset of flutter that oscillations will begin to build up 
prior to the blade reaching a critical azimuth position, then decay as the blade 
moves beyond this point. This build up and decay means that the vorticity shed 
due to the oscillations will be contained within a double azimuth region on either 
side of the critical azimuth position. Assuming this region to be small allows the 
wake system to be two-dimensional. 

The lift deficiency function resulting from the two-dimensional wake model is 
compared with earlier results obtained for a helicopter in hover and fixed wings. In 
the limidng case when the advance ratio is very small, Shipman and Wood’s lift 
deficiency reduces to Loewy’s lift deficiency function. Also in the limit as the 
wake spacing is very large, this function reduces to Theodorsen’s lift deficiency 
function . 



' Shipman, K. W , and E. R. Wood " A Two-Dimensional Theory for Rotor Blade Flutter 
in Forward Flight." Journal of Aircraft, Vol. No. 8, Number 12, pp. 1008-1015. 

December 1971 
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The theory was applied to bending torsion flutter for the tip segment of a rotor 
blade. The results followed the normal trends of having the flutter velocity decrease 
as the blade center of gravity was moved aft and of having the flutter velocity 
increase as the blade stiffness increased . 

The previously shed wakes are destabilizing so that they reduce the flutter 
velocity. The buildup and decay produces wakes that are essentially centers of vor- 
ticity Thus their position with respect the reference blade will be sensitive to the 
advance ratio which determines the vertical spacing. It was found that as the 
wakes were brought closer by reducing the advance ratio and/or inflow ratio, the 
flutter speed is decreased [29]. 

The phenomenon of airfoil flutter is caused by energy extraction out of the 
wind stream by means of airfoils which are free to oscillate in pitch and plunge. 
As described by Duncan [19], this effect can be easily demonstrated in a wind tun- 
nel if a device is built which allows to vary the phase relationship between pitch 
and plunge. Phase angles near 90 degrees produce maximum energy extraction so 
that the device operates as a flutter engine. McKinney and DeLaurier [20] demon- 
strated that this effect can be used to convert wind energy into mechanical energy. 
Scherer, McKinney and Delaurier used Theodorsen's flat-plate theory for perfor- 
mance estimates of the flapping foil propulsor and of the wingmill, respectively. 

It was shown by W. Schmidt, that optimum energy extraction occurs for har- 
monically flapping airfoils upstream of a stationary airfoil. There is a need for 
aerodynamic analysis methods which enable systematic studies of the effect of air- 
foil geometry, incidence angle variation, amplitude of oscillation and airfoil interac- 
tion on the aerodynamic forces due to incompressible flow over airfoils or airfoil 
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combinations executing pitch or plunge motions. A panel method was developed in 
order to do parametric studies of these effects. In the past, a number of investiga- 
tors have solved the steady tlow problem using source and vortex paneling, the 
most prominent ones being Hess and Smith [21] A few authors have extended this 
approach to the case of unsteady motion of single airfoils, notably Basu and Han- 
cock [22] and Kim and Mook [23], 

Basu and Hancock [22] presented a numerical method for the calculation of 
the pressure distribution, force, and moments on a two-dimensional airfoil undergo- 
ing an arbitrary unsteady motion in an inviscid incompressible tlow. The method 
was applied to (i) a sudden change in airfoil incidence, (ii) an airfoil oscillating in 
high frequency and (iii) an airfoil passing through a sharp-edged gust. 

Teng [24], developed a computer code for the numerical solution of unsteady 
inviscid incompressible flow over an airfoil. Pang [25] extended the work of Teng 
and developed a computer code to study unsteady airfoil interference effects. The 
technique uses the well known panel methods for steady flow and extends it to 
unsteady flow by accounting for the continuous shedding of vortices into the trailing 
wake. 



29 



III. DRAG REDUCTION DUE TO 



BLADE FLAPPING 
( KATZMAYR EFFECT ) 

A. INTRODUCTION 

In this chapter, a discussion will first be presented of the physics of dapping 
airfoils which explains the bird's ability to generate a forward thrust by means of 
wing flapping. This will be followed by a description of two methods capable of 
predicting the pressure distributions, lift, pitching moment, and drag or thrust on airfoils 
which execute harmonic oscillations in plunge or pitch. 

The first method is based on the generalization of the well-known steady-state 
panel method to the case of unsteady airfoil motion. It is based on the assumption 
of inviscid, incompressible flow. The second method is based on the classical The- 
odorsen theory for oscillating flat-plate airfoils in incompressible inviscid flow. The 
unsteady panel method permits the investigation of finite-thickness airfoils of arbitrary 
geometry which may execute rather general unsteady motions ( ramp and oscillatory ). 
In contrast, Theodorsen’s analysis is restricted to thin airfoils (flat-plates) and to small 
amplitude oscillations. The chapter ends with a comparison of results computed by both 
methods. 
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B. THE PHYSICS OF FLAPPING AIRFOILS IN LOW-SPEED FLIGHT 



The phenomenon of thrust generation due to wing flapping or due to steady flight 
in a sinusoidal gust can be explained in a rather elementary way. Consider the latter 
case of flight through a sinusoidal gust, as shown in Figure (3.1). The relative wind 
vector changes direction between a maximum positive and negative angle of attack. The 
resultant aerodynamic lift is perpendicular to the instantaneous wind vector. If we as- 
sume the viscous drag to be small, then it is easily seen that the sinusoidal gust gen- 
erates a sinusoidally varying thrust. 

A similar effect is generated if the airfoil executes a sinusoidally varying plunge 
motion ( wing flapping ) about an otherwise steady flight conditions. This explanation 
of thrust generation due to wing flapping or due to flight through a sinusoidal gust was 
first advanced by Knoller [10] and Betz [11], Its experimental demonstration was 
accomplished by Katzmavr [ 7 ] and is therefore generally referred to as Katzmayr 
effect 1 

Since it is well known that every change in airfoil incidence causes the shedding 
of a starting vortex, a more precise explanation of the Katzmayr effect has to take into 
account the continuous vortex shedding due to wing flapping or due to flight through 
a sinusoidal gust. 



1 Platzcr. M. F., Neacc. K. S , and Pang. C. K.," Aerodynamic Analysis of Flapping Wing Propulsion", 
page 1, AIAA paper No. 93-0484, 31st Aerospace Sciences meeting. Reno. NV. January 11-14. 1993. 
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Figure (3.1) Propulsive force on plunging airfoil 
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Similarly, the quantitative prediction of the induced oscillatory forces ( lift and 



thrust ) had to await the development of a complete unsteady aerodynamic theory. 
Such a theory was first accomplished in the 1930’ s by Theodorsen [3] in the United 
States and by Kiissner [16] in Germany. It can be found in the textbooks on aeroelas- 
ticity, for example the book by Bisplinghoff. Ashley, Halfman [41], Therefore only 
the final formulas will be given here and the reader is referred to those texts for 
additional details. More recently, a second approach to the determination of the 
unsteady aerodynamic forces on oscillating airfoils became feasible with the advent of 
high speed computers and the development of efficient numerical solution methods. The 
fundamental building blocks of this method are presented in the next two sections. 

C. PANEL CODE FOR STEADY INCOMPRESSIBLE INVISCID FLOW 

In potential flow theory, the flow field around an airfoil may be represented by 
the velocity potential. Considering contributions from the free stream flow and the 
source and vorticity distribution, the total potential may be constructed: 

d) — d^p ■+■ “i" Oy (3.1) 

where 



d>p = U [ x cos ( a ) + y sin ( a ) ] 
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The source distributions q(s) vary front panel to panel, while the vorticity strength 
Y (s) is assumed constant for all panels. The value of representing the flow past an 
airfoil by surface singularity distributions lies in the fact that these singularity distri- 
butions automatically satisfy Laplace's equation, the governing flow equation for invis- 
cid incompressible flow: 



0+0 = 0 (3.3) 

\\ yv 

Since Laplace's equation is a linear homogeneous second order partial differential 
equation, the superposition principle used in Equation (3.1) holds. The boundary con- 
ditions include flow tangency at control points (mid-points of panels) and the Kutta 
condition at the trailing edge, requiring equal tangential velocities for the first and last 
panels. By evaluating the integrals along the airfoil surface, the potential may be de- 
termined at any point in the flow field. Each point is defined at a radius (r) and angle 
(0) from a chosen reference point on the airfoil. 

The reference point in this study is the leading edge. The airfoil is represented 
by a number of defined points, called nodes [34], More points produce greater resolution 
and better accuracy. One hundred to two hundred points are usually sufficient, with 
the larger number used for more complicated airfoil shapes or more involved calcula- 
tions. The lines connecting these nodes are the panels. There are n+1 nodes with first 
and last node overlapping. 
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Figure (3.2) depicts die panel geometry Numbering starts at the trailing edge, then 
resses along the lower surface, leading edge, and upper surface, and ends at the 
ng edge. The unit normal vectors (n,) are perpendicular to the panels and directed 
aid from the airfoil surface The unit tangential vectors ( t, ) are parallel to the 
Is and the positive direction is defined with incrcasiim numbering (n to uH ) 
panels may vary in length, with the exception of the tlrst and last panels, which 
be equal in order to use the Kutta condition at the trailing edge. 
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Figure 3.2 Panel Method Geometry 
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The required input consists of the number ot nodes on the airfoii surface, the 
coordinates of the nodes referenced from the leading edtie. and the angie of attack in 
degrees. The program produces normalized velocities and pressure coefficients at each 
control point as output. The use of influence coefficients leads to a straightforward 
procedure for programming the equations 

An aerodynamic influence coefficient is defined as the veiocitv induced at a field 
point by a unit strength singularity distribution on one panel. For the two-dimensional 
steady How problem, the following influence coefficients are needed 1341 

A n n normal velocity comDonent induced at the i 1 panel control point by unit 
strength source distribution on the j th panel 



A n = sin (0 - 0 ) In 
'J 27i. 1 •' i 





(3.4) 



r 



1/2 



A l jj tangential velocity component induced at the i 11 panel control point by unit 
strength source distribution on the i' h panel 



A.. = 



ij 




(3.5) 



ii 



= 0 



B n lt ; normal velocity component induced at the r !l panel control point by unit 
strength vorticity distribution on the j' !l panel 



1 r | 

B n = cos(0 -0)ln " -sin(0 -0)(3 
ii 2 ^' i j r i .r u 



(3 6) 



= 0 



i =J 
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B l n tangential velocity component induced at the r 1 panel control point by unit 



strength vorticitv distribution on the i :i panel 



B f - cos (0 - 0 ) 13 - sin ( 0 - 0 ) I n 1 
•J 2 k ' ' ’i 1 ' r 



l = i 



(3 7) 



u 



- 1/2 



where the geometrical quantities, depicted in Figure (3 3) are defined bv 



r . = (xni - \ )“ - (vm - v ) . \m 

ii v t i ^ i - i i 



x + x 

i i-i 



vm = 



V T V 
- 1 



0. = arc tan 



I r 1 



\ , - X 



S) 



( ym - y. . ) f ym -y3 

n J i + 1 x h 7 i J j 

p = arc tan i — l-arctani 

'J l xm - x , J Ixm. - x J 

v i j + l y v i y 

The first boundary condition requires How tangencv at the control points 

(V n ), = 0 , i = 1, 2 n (3 9) 

In terms of inlluence coefficients (with V = 1 ), 



I 

j = i 



\ n 

A..q. 
ij .1 



+ v y B 

z- ij 
j = i 



sin (a - 0 ) = 0 



= 1 ? 



, n 



(3 10) 



The second boundary condition is the Kutta condition, which states that the pressures 
on the lower and upper panels at the trailing edge must be equal if the flow is to leave 
the trailing edge smoothly. 
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I'igiirc (3.3) Relationship Between Geometrical Quantities 
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Using Bernoulli's equation, the pressure coelllcient can be expressed as 



I - 






The pressure equilibrium also implies equal velocities lor incompressible How Since 
thenounal velocities are taken to be zero, the boundary condition may now be staled as 



wlicie the negative sign is strictly due to the adopted convention of positive tangential 

velocities in the direction of increasing node numbering Since the How is positive to 

the light, the panel downstream of the front stagnation point will have negative values 

for computational purposes only. It is important to note that not all the lower surface 

panels have a reversed sign, only those downstream from the stagnation point. This is 

especially significant for non- symmetrical airfoils or any airfoil at an angle of attack. 

In terms of influence coefficients, the normalized equation becomes . 
n n n n 



liquations 3.10 and 3 13 represent a linear algebraic system of (n+| ) equation and 
(nt I) unknowns. The unknowns are the source strengths which vary from panel to 



(V), - (V*)„ 




panel ( q ( 



q n ) and the vorticitv strength y. 
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Expanding and rearranging Equation 5 ! 0 for an example airr'oi! of n = 7 3 noaes and 
panels results in 



"•Aii 


4 

r l 
r 1 

1 




- B : ,A - .. 


. - B'. 1 - 

I . A 


-'in (a -f .) 


-A 


- A; : q : - .. 


-WBt, 




• " B :. = 


'in tot - d, ) 



A-, ", q | - A2, - ... - y ( BIT , B'2, , - . - B1‘, = ->tn (a - A-,, 

The equations now readily lend themselves to solution in matrix form Recasting 
with simpler notation, the A''\, terms (coefficients of q, ) mav be renamed a., ana the 
sum of ail B n „ terms in parentheses (coefficients of y > renamed a, nri . where i = 1 . 

2. .. . n and i = 1 . 2. n. The terms on the right side of the equation may 

be renamed b, . The f n-+- 1 ) equation, or in this example, the 7 -T 11 equation comes 
from Equation 3.13 in a similar manner: 

( a; , t At. , ) q , - ( A 1 , . - At. ,) q,. ' (A‘ . At. q,, 

-7 I ( B* | | - B 1 -,, ) - ( B' j . t- - ( B‘ B 1 -,-. )] 

= cos (a-Gj) - cos ( a - 673 ) (3.15) 

The coefficient q[ may be renamed a-r 4 All of the B l terms in the brackets together 
form the coefficient of 7, now renamed a 74.74. The entire right side of the equation 
constitutes the new term b-4 
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Finally expressing this system can be written in concise matrix form 



a 1 . n + I C 1 1 

a :.n-i V 



b. 



a , a 
n. 1 a. 



a , . a 

n + 1 . 1 



n. n - i q, 



n + I , n + 1 L ! _ 



b 

■I 




(3 16) 



This system is solved in the program using a Gaussian Elimination subroutine. With 
the values of the q, and y known, the velocity at each panel control point may be 
calculated 

n 

V. = ^ Ap^i -cos(a-9 t ) , i = 1, 2, n. (3.17) 

J = I 

The total velocity is equal to the tangential velocity due to taking the normal velocity 
to be zero. 
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D. UNSTEADY PANEL CODE 



The extension of this method to an airfoil experiencing .. rapid change in 
angle o» attack can be accomplished as follows It is important to iccaii Pranatl’s 
eariv tlow \isuaiization experiments which demonstrated the generation of a 
starting vortex off the trading-edge whenever the angle of attack was changed. 
Hence continuous changes in angle of attack produce a continuous snedditm of 
vorticitv into the trailing wake. This shedding can be explained bv the Heimiioltz 
cortex theorem wiiicit requires that anv change in the circulation around the airfoil 
must be matched bv the appearance of an equal counter\orte\ at the trading 
edge in order to achieve constancv of the total circulation .in tne flow field. 

The vortices shed from the trailing edge move with the fluid particles of the 
surrounding fluid and hence are swept downstream with a speed essentially equal 
to the free-stream speed. Therefore these vortices will stav close enough to the 
airfoil lor a finite time to influence its pressure distribution. Any unsteady airfoil 
theory therefore must describe this vortex shedding process. It is this feature 
which distinguishes unsteady airfoil theory from its steady counterpane The 
above considerations suggest that the airfoil can be modeled by similar source 
and vonex distributions as in the steady case. 



* Platzer. M. F. Class Lecture Notes. Navai Postgraduate .School. Monterey. 
California. Sep. p)9.' 



42 



Hence source panels are again placed on the airtbii surface together with a 
single vortex distribution which is the same for each panel An arbitrary change in 
angle of attack then is subdivided into small step changes such that a starting 
vortex is shed into the wake each time a step change occurs and the flow tangencv 
and Kutta conditions are enforced at each time step. 

In addition to the n airfoil panels a wake panel is assumed to be attached to 
the trailing edge such that, after each time step, the vorticity of the wake oanel is 
assumed to be concentrated into a single point vortex which detaches from the 
airfoil and starts moving downstream with the tree-stream veiocitv Vs one wake 
panel detaches from the airfoil a new one is created and a wake of point vortices 
is created As in the steady problem n unknown source strengths and one unknown 
vorticity strength on the airfoil are introduced. 

In addition, there is now the unknown vorticity strength of the wake panel 
whose length and orientation are also unknown. Hence there are three additional 
unknowns which need to be determined by a proper set of equations; (a) the 
flow tangencv conditions provides n equations, (b) the Kutta condition provides an 
additional equation, (c) the Helmholtz theorem provides an equation for the vortic- 
ity strength on the wake panel. 

Therefore two more equations are needed to determine the length and orienta- 
tion of the wake panel. These are obtained by making the following assumptions, i.e.. 
(d) the wake panel is oriented in the direction of the local resultant velocity at its 
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md-poim. tei the letigtn of the '.\ake panel is propomonai to the magnitude of the 
local resultant \eiocity at us mid-point and to the time-step The Kutta condition 
and the last two conditions are nonlinear equations and therefore necessitate an 
iterative solution procedure Figure (3 4). from Ref 1251. shows the panel method 
representation for unsteady flow 



44 



i (elmlioltz s theorem 



Vortex Shedding nt lime Step t^ 



f k 



k-l 



Panel j 

Source Distribution (q.) k 
Vorticity Distribution y. 




( V A 



lanO. 

k 



< V A 

C'wV 




'V “ Ck~'i.|) V((U„) k ’ k- (VJ„ ! | 



Figure (3.4) Pniicl Methods Representation for Unsteady Flow 
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Z. THEODORSE.VS ANALYSIS 

As stated before. Iheoaorsen's theory for osculating riat-piate airfoiis in iiimsckI 
ncomnressible f:o\v is well documented. Therefore, e summarize .'ere tniv the 
final formulas. The lift on an airfoil oscillating at freciuencv co is given bv 

L - - per V ; Li, (h. v b) - \ L. t - < ! 2 - i i Li, , a ; i 2.18) 

\ nere 

L h = ! - , 2 i C(k) ) k . 

L a = 1/2 - ( i / k ) [ I - 2 C(k) ] - ( 2 / k ; ) C (k). 
k = co b / U ( b is the semi-chord and "a" is the elastic axis ) 

C(k) is the lift deficiency function and is given by 

H, (2) (k) 

C(k) = — = F (k) - i G (k) (3.19) 

II, (k) + i (k) 

where H“ n (k) is the Hankel function of second kind and is given in terms 
of Bessei functions as H“ n (k) = J n (k) - \' n (k ) 
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F. GARRICK’S PROPULSIVE FORCE ANALYSIS 



Garrick [5] applied Theodorsen's theory to determine the thrust and drag 
generated by an airfoil which executes a sinusoidal motion in pitch or in piunge. 
He found the following formula for the average propulsive force 



P x = rcpbco" a, h~ - (a, - b,) a ( - 2 (a 4 - b ( ) a h 



(3.20) 



where a. F - - G" 



a 2 - b - ! a i [ I k - - 1 - ( 1/2 -a) - J -1/4 - (1/2 - a ) F - G ' k ! 



a 4 = b ! a j [ - sin ( <p 2 - <p n ) + (1/2 - a ) cos <cp 2 - <P 0 ) ] 

- (F/2) cos ( cp 2 -<p 0 ) + (G/2) sin ( <p 2 - cp 0 ) | 
b 2 = b 2 [ - a/2 - F/k 2 + (1/2 - a ) (G/k) ] 



b 4 = (b/2) [ (1/2 G/k ) cos ( (p 2 - tp 0 ) + (F/k) sin < <p 2 - (P 0 )] 

For pure plunge, the average propulsive force is reduced to 

P N = TtpbcoV F 2 (k) +G 2 (k) i (3.21) 



For pure pitch, the expression for the average propulsive force is given by : 

3 2 ^ 

P,x = Ttpb CO a n [a 2 + b 2 ] (3 22) 



where a 2 and b 2 are defined in Equation (3.20). 
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The expression for the average drag in Equation (3 20) vt Garrick's reDort 
'o7 [5], was found to be m error The correct expression as given in Garrick's 
subsequent review [32] is 

^ 'll i _ 1 r ! G i 

D = VC - • - \ - -F -x- - -x . 1 3.23) 

2 2 • / 2 2 v - 2 • l; . 

k x 

w here 

C" - F" - Cr and \ - ~ . p. o". 'ir l v ' 

Cm. comparison of results 

in a recent thesis. Riester [43] presentea a uetaiied comparison with 
Theodorsen s predictions of the lift and pitching moment computed with the un- 
steady panel code. Generally, good agreement was found, thus lending credence 
to the reliability of the panel code. In the present work, we are interested in 
the prediction of the thrust and drag forces generated by an oscillating airfoil 
Calculations were performed to ascertain the panel code's capability . Figures 
(3.5) to (3.7) show the comparison between Garrick's predictions and the nanel 
code results. These figures present plots of the average propulsive force coef- 
ficient vs the non-dimensional plunging amplitude of oscillation. HBAR (HBAR 
= h/c ), for three reduced frequencies. k = 0. 1 . 1 and 2. The two methods are 
in good agreement for small plunge amplitudes. As might be expected, the 
agreement starts to deteriorate with increasing plunge amplitude. The panei code 
was applied to an MAC A 0009 airfoil using 100 panels. 
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Figure (3.5) Comparison of Panel Code Resulls ,viH, Garrick ( k = 0.l ) 
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Figure (3.6) Comparison of Panel Code Results with Garrick ( k= I ) 
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f igure (3.7) Comparison of Panel Code Results with Garrick ( k = 2 ) 
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H. ESTIMATED POWER REDUCTION DUE TO UHC FOR OH-6A 

HELICOPTER 

The ultimate objective of this research was to cietermine tne effect or unsteadv 
aerodynamics on the measured power values for tne OH-6A heiicoprer. Ve are 
especially interested in knowing whether these effects couid produce the power 
changes recorded during the NASA - Army open loon flight tests in the lligitt 
regime from hover to 100 knots. 

In order to estimate the change in torque due to propulsive force ( foot- 
pounds per blade) due to HHC. the blade span was divided into sections As- 
suming hover conditions, the reduced frequency at each blade section was cal- 
culated The panel code was used to compute the propulsive force coefficient. 
Cd, at each blade section assuming pure plunge motion with a plunging am- 
plitude of 1". The blade deflections, obtained from Wood et al [27], were then 
used to estimate the real value of the plunging amplitudes at each blade segment. 

Finally, the net blade torque could be calculated by considering the propulsive 
force at each blade section and multiplying this value times the rotor radius, then 
summing over the length of the blade. Table 3.1 . from Ref. [42], shows the 
OH-6A helicopter parameters that were used to estimate the propulsive forces. 
The most important parameters are: number of blades. Q =4. rotor radius (ft), 
R= 13.2. blade chord (ft), c = 2b = 0.57. rotor speed (rad/sec). Q =49.2. and 
helicopter weight (lb), G = 2550. 
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TABLE 3.1 THE OH-6A HELICOPTER ROTOR PARAMETERS : 



Rotor radius (ft) 
Number of blades 
Rotor speed (rad/seci 
Blade chord (ft) 



4 

50.26 

0.57 
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Table (3 2) shows, at eacti blade section r R. the reauced frequency. k(r). 
local speed. LTr), the computed average propulsive rcrce. P. wri ana the torque. 
T, n w> due to the propulsive force at each biaae section, summing me esti- 
nated toraue at each blade section, the resultant toraue uue to the pronuisive force 
per blade was found to be 8 11 ft-lb assuming .. pure munging motion with 

constant piunging amplitude of 1 " u each blade - ecu on. 'Food et ai [27], ap- 
plied Garrick's equation and considered the effect of the second to twelfth har- 
monic loads He found, for 2P harmonic < HHC n ’ ? 30 degree rnase at 60 

knots) that the power gained is about 11.3 horsepower “or more uetaii. the 
reader ts referred to tabled of Ref. [271 Figures 7 > to • :! now me oiade 
deflections for several harmonics used by this reference mown m figure (3 12) 
is the average drag coefficient. Cdav. \ersus reauceu frequencv as obtained with 
the panei code for a plunging amplitude of 1" The airfoil used was VAC A 0015 
(OH-6A airfoil) The propulsive force. P^ , is given by 

P xh = 0.5. p . U 2 C dr. C d (3 24) 

and the torque. T\h. at each blade section is given bv 

T \h = p \h r 1 

where U is the local speed at each section . C d is the propulsive force coef- 
ficient as obtained by the unsteady panel code and r is the radial distance of 
each blade section For the OH-6A. the angular freauencv . Cl = 50.26 rad/sec. 
therefore the power per blade is found to be 0 7- 1 horsepower as giv en bv 

P = Txh Q . (3.26) 
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TABLE 3.2 LOCAL REDICED FREQUENCY. K p (r). PROPLLSIYE 
FORCE. P %h ( r) AND rORQl'F.. T xh (r) . 
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k (r/R) 
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Figure (3.10) Blade Deflection vs Radial Station at 4P 
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Figure (3.11) Blade Deflection vs Radial Station at 5P 
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versus Reduced Frequecy ( for 1 " plun^ii)” amplitude ) 




Figure (3.12) Cilav vs kp ( Plunging Amplitude 1" ) 
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IV. ANALYSIS OF WAKE 



INTERFERENCE EFFECTS 



A. INTRODUCTION 

The results of Chapter III show that propulsive forces, generated soleh 
by plunge oscillations of an airfoil are insufficient to expiain the power 
reduction observed and measured on the OH-6A helicopter Therefore, in this 
chapter, we examine whether the power reduction might be caused by favorable 
wake interference effects due to the wakes shed from preceding blades. 

To this end, the extension of the unsteady panel code to the case of two 
airfoils [25] was adopted and applied to the study of the wake interference effects 
between two helicopter blades. Wake interference effects were studied in the 
1950’s by Loewv [4] by extending Theodorsen's analysis. Loewv treats an 
infinite number of layers of shed vorticity placed beneath the rotor at a given 
wake spacing and developes a modified lift deficiency function His interest 
was in the possibility of wake-induced blade flutter, which restricted his analysis 
to the determination of lift and pitching moment changes due to wakes shed 
from preceding blades. A closer inspection of Loewv's paper suggests that it 
should be possible to use the same approach to examine drag and propulsive 
force changes on the reference blade due to wakes shed from neighbounng blades. 
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Further, it is desirable to earn.' this one step further and investigate the 
ease of a single wake shed from the preceding blade which interferes with the 
wake shed from the reference blade 1 This ease allows a direct comparison 
with the unsteady two-foil panel code. 

Also, there are regimes of helicopter flight where only the first iaver of 
shed vorticitv is of consequence. In the following sections we first describe the 
two-foil panel code and Couch's modification of the Loewv anaivsis This is 
followed by a presentation and comparison of the maior resuits. 

B. UNSTEADY TWO-FOIL PANEL CODE 

The extension of the single airfoil code ( previously described ) to the case 
of two airfoils requires no new building blocks. However, the two-foil analysis 
requires the introduction of five frames of reference, namely two moving local 
frames of reference which are attached to the two airfoils, two frozen local frames 
of reference, and the inertial frame of reference. Furthermore, it requires the 
satisfaction of the two Kutta conditions which are coupled non-linearlv The 
solution requires an iterative procedure to compute the two vorticitv distributions. 



1 A detailed development of this theory is derived and presented 
in a recent NFS thesis bv Mark Couch. Ref. [40 1. 
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It also necessitates the creation of a subroutine which transforms ail coor- 
dinates in either of two local frames of reference to the global frame of reference 
Finallv. it requires the extension of the influence coefficient concern to include 
the effect of the second airfoil with its own wake anu it reciuires the introduction 
of an additional influence coefficient, namelv that on the wake element due to 
the wake element from the other airfoil If each airfoil is modelled with X 
panels, then this produces a 2X T \ ( X-3 ) matrix wnich .s HiDseuuentiy -oived 
by Gauss elimination 

Pang [25] was only able to \erify this code w .ompanng me computed 
pressure distributions with previous work by Giesing [531 Hence it was nec- 
essary to further evaluate the code by comparing its output with otner known 
solutions before applying it to the oscillatory blade interference problem. These 
comparisons are given in section C. 

The main parameters used in this code are the airfoil type, number of 
panels. N. the relative locations of the two airfoils, initial angle of attack, am- 
plitudes of oscillation (pure pitch, pure piunge or both), frequency of oscillation 
for each airfoil, the pivot points (for pitch cases), rise timet for ramp motions ). 
These input parameters are in the file "fort 1" 

This program produces an extremely large amount of output to the screen. It 
is convenient to write the output to a file on a Unix based machine during 
program runs. The output file can get very large for long program runs (when 
the time history is long). 
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For this reason, the logical variable 'output" was added to the input file m two 
inodes true" or 'false" When this variable is set to false", the screen output 
is not printed and the required outputs are redirected to output files 

The following list describes the output riles and the data thev contain 
ton 2 This is for user supplied airfoil coordinates, if desired 
fon 3 This rile contains the global coordinates of \F1 at each time step 
fon. 4 This rile contains the global coordinates of AF2 at each time step, 
fon. 7 This rile contains the lift, moment, and drag coefficients for Doth airfoils 
fon. 8 This rile contains the pressure coeff. for AF1 at each time step, 
fon. 9 This rile contains the pressure coeff for AF2 at each time step . 
fort. 10: This rile contains the first airfoil’s core vortex (wake) positions, 
fort. 1 1 This tile contains the second airfoil’s core vortex (wake) positions 
fort. 12: This file contains the angle of attack at each time step, 
fort. 13. This file contains the y-axis translational motion at each time step. 

C. LOEWY’S ANALYSIS 

Loewv [4] developed a two dimensional model for the unsteady aerodynamics 
of the blades of a hovering rotor, including the effect of the shed wakes. Figure 
(4 1). from reference [40], shows his two-dimensional model of the helical 
wakes. The wake in the case of a multi-bladed rotor in hover consists of 
helical vortex sheets below the disk, one from each blade. 
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Assuming sinusoidal motion, the downwasn over me reference airfoil. V, and 
its vorticitv (either bound or shed) may be written, respectively as 

V = v\e'" U . v = v e " l (4 1 

•a ■ a c 

where the barred quantities are complex The induced velocity at a point ; on 
the airfoii resulting from an element of vorticitv of strength. y l(J at a general 
point in the wake is given by [4] 



where Q is the total number of blades and h is the vertical distance between 
successive rows of vorticitv. n is the number of revolution index, and q is the 
number of blades index as shown in Figured 1 The wake spacing, h is defined 




4 2) 



dw < x.n 



2 — J ( x — 4 > - ( nQ - q) (h) " ) 



as [4] 



h 




(4 3) 



where u is the inflow velocity ( time-averaued normal to the disk) 
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I'igitrc (4.1) Aerodynamic Unsteady Model for a Mnlti-Bladed Rotor. 



The total downwash over the airfoil can be expressed conveniendy by summing 
the integrals involving the bound voracity and the row of vorticitv in the plane 
of the reference blade with those for the rows of vortices below the plane of the 
rotor 
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The first two terms on the right-hand side of this integral equation are those 
that arise in the classical unsteady airfoil theory’ of Theoaorsen | 3 I. the third 
and fourth terms contribute the influence of all the voracity beiow the oiane 
of the rotor disc Loewv solved the integral equation using Sohngen's inversion 
formula [31], and obtained his lift deficiency function. C tk.m.ii). wnich is given 
bv 

C (k.m.h) = F (k.m.h) - i G (k.m.h) (4 5) 

where k. is the reduced frequency (detinea bv o l ' ). m is given hv he 

ratio o/D. h is the wake spacing defined before, and the real ana imaginary 
pans of Loewy's lift deficiency function. F . G are given by [4] 



J.ctA- (Y. -pJ,)B 



F’ (k.m.h) = 



A' -r B" 



(V, -r pj. ) A -T J aB 



G' (k.m.h) = 



4" - B" 



w here 



A = J a + Y -J p 

1 O O 



-B = Y. -r J,P- J a 
1 r o 



and 



kh -kh 
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kh ... 

e - 2 cos (2^m) -e 
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Here it can be seen that as the spacing between rows of vorticitv, h, 
becomes infinite. P approaches l : and Ci’ approaches G. Also as k tends to 
infrnity, C’(k,m,h) approaches C(k) When h is zero, all the vorticitv lies in 
the plane of the airfoil, singularities result, and the meaning of anv obtained 
results, physically speaking, is not clear [4 J Loewv [4 J, and Couch [40] showed 
plots for F and Ci versus reduced frequency k. Loewv treats an infinite 
number of wakes and Couch explores a finite number including the special case 
of a single wake. These plots show the effects due to wake phasing, m , and 
wake spacing, h A case of special interest is the wake shedding from one 
preceding blade only, for this ease l.ocwv's lift deficiency function, C‘ . reduces 
to the following lift function. C*, obtained by Wood and Couch [ 40] 



where 



C* 

F* 



= F* - i G * 
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The propulsive force coefficient can be calcuiated and the results, obtained by 
Wood and Couch [40], using F* and G* are : 



For pure plunge 
For pure pitch 
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D. EVALUATION OF THE TWO-FOIL PANEL CODE 



Several analyses are available in the aerodynamic literature which iterant a 
: Tined evaluation of the two-foil panel code These are: 

’round Effect on Airfoil Lift in Steady Flight 

S. . .. ika et al [44] analyzed the around effect on a 2-D flat elate using 
the method of conformal transformations Figure 4 2a shows the panel code 
results, Figure4 2b reproduces Tomotika results. It seen that the two metnods 
are in good agreement. Figured 3 a and Figure 43b show the comparison be- 
tween the result obtained by Weglev [46] for a i 0 0 o thick RAE101 airfoil in 
ground effect and the panel code prediction for a NACA 0010 at zero angle of 
attack. The two results are again seen to be in good agreement considering the 
fact that the RAE101 airfoil (whose coordinates are not known ) was approxi- 
mated by a NACA 0010 airfoil. 

2. Bip^-.e Effect 

Glauert [45] defined a factor B which gives the lift reduction on the airfoil 
in an unstaggered biplane configuration compared to the lift on a monoplane at 
the same incidence. Table 4,1 gives the correction factor for several non-dimen- 
sional biplane spacings. It is seen from this table and from Figure (4.4) that the 
two results are again in reasonable agreement. 
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CL/CIo 



C/CIo vs h/c 




h/c 



Figure 4.2 a CI/CIo vs h/c ( Panel ( ode) 




Figure (4.2 b) 0/0 o vs h/c ( From Ref. [44| ) 
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I'igurc (4. 3;i) Local Velocity Distribution, V„/V vs X/C (Panel Code) 




r inure (4.3 b) Local Velocity Distribution, V,/V vs X/C (after Wegiev (46|) 
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0 

Figure (4.4) Correction Fnctor, B vs (he Clioni-gap Ratio h/c 
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TABU*: 4.1 CORRECTION FACTOR FOR A BIPLANE 



h/c 



B ( Code ) B ( ref. |45| ) 



0.50 


0.770 


1 

0.730 


0.57 


0.847 


0.800 


1. 00 


0.902 


0.855 


1.25 


0.942 


0.895 


1.50 

_ __ 1 


0.97 


0.920 

J 



Another ground interference study is due to Bagley [46], As can be seen from 
the comparison of the panel code (Figure 4 5 a) with his predictions ( Figure 4.5 
b ), in ground effect, the lift is reduced over the upper surface and increased 
over the lower surface As the airfoil approaches the ground, the flow between 
the airfoil and ground is reduced and hence the pressure is increased. 
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Cp vs X/C (SS.il/c- 5. AOA=l I cleg.) 




X/C 



Figure (4.5 a) Pressure Distribution vs X/C ( Panel Code ) 




'•Sme (4.5 b) Pressure Distribution vs X/C ( Daglcy) 
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E. ANALYSIS OF OSCILLATORY BLADE INTERFERENCE EFFECTS 



Having determined in Chapter ill that pure flapping of a singie heiicomer 
rotor blade is unlikely to produce a drag reduction sufficient to explain the 
OH-6A HHC power reductions [ I ]. soielv. due to me Katzmayr effect, we 
investigate in this section the possibility of favorable blade/wake interference 
effects Consider the arrangement of the two airfoils shown in Figure (4 6) The 
wake shed from airfoil -2 impinges airfoil ~ I. if the \erticai spacing between 
the two airfoiis falls below a certain \aiue. depending on the antoiituue of 
oscillation of airfoii #2 This case cannot be analyzed with the unsteauv panel 
code because the cortices start to penetrate the second airfoil or oecause the 
wakes shed from the two airfoils come in contact with each other, as shown 
in Figure (4 7). 

On the other hand, if the vertical non-dimensional spacing is two or more 
no problems occur A second important parameter is the horizontal spacing be- 
tween the two airfoils. This parameter controls the phasing between the two 
wakes. Hence the unsteady two-foil code permits the analysis of Loewv's ar- 
rangement for the special case of two wakes. This case is considered by Wood 
and Couch [40] in closed form and presented in Ref [39] Furthermore, the 
single airfoil case is recovered if the vertical spacing is chosen to be large 
because the effect of the first airfoil on the second must diminish with increas- 
ing distance between the two foils. 
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rignre (4.6) Wnke Pattern (Airfoil #2 Plunging at High reduced Frequency 
and Close (o Airfoil # I ) 
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Figure (4.7) Wake Interference Between Two Plunging Airfoils in Close 
Proximity 
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Figure (4 8). and Figure (4 9) show the computed lift and drag (thrust) co- 
efficient for airfoil n 1 in the presence of airfoil ~ 2 Both airfoils are oscillating 
in plunge with an amplitude of 0 14 and a reduced frequency' of 0 0617 The 
vertical spacing between the two airfoils is 200 In this case the lift and thrust 
coefficients are found to be identical to the ones obtained from the single airfoil 
code, thus showing that a distance of 200 is sufficient to recover single airfoil 
results 

In Figures (4 10) to (4 21). on the other hand, the vertical spacing is reduced 
to 2 For these calculations, the time step is 3 to .5 time units The piunging 
amplitude is again 0 14. the reduced frequency is 0 0617 The twelve figures 
4 10 to 4 21 show four different cases of phasing between the two wakes. In 
Figures (4 10) to (4 12), the wake pattern, lift and drag are shown for zero 
phasing, m = 0 ., Figures (4. 13) to (4.15) show a phase of m = 0.20833, Figures 
(4 16) to (4.18) show the case of m = 0.25, and Figures (4 19) to (4 21) show 
the case of m = 0.5 . 

A closer inspection of these figures reveals the strong effect of the wake 
phasing. For m = 0, airfoil # 1 experiences a reduction in thrust compared to the 
single foil, whereas for m = 0.5 airfoil #1 experiences a significant increase in 
thrust. These results are summarized in Figure (4.22) where the single airfoil 
results are also shown as well as the predictions by Wood and Couch [40] 
when specialized to the case of two interfering wakes only. The agreement be- 
tween the two predictions is encouraging. 
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Figure (4.8) 



Cl vs T* 



( plunge, m = 0, It = 200 ) 
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Non-Dimensional Time iT*j 




Figure (4.9) 



Ctl vs T* ( plunge, 



m = 0, h = 200 ) 
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igurc (4.10) 



Wake Pattern ( plunge, in = 0, li = 2 
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Gl,C12 vs T' (kl=k2=.0617,m=0,h=2,hbar=. 14) 
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Figure (4.11) Cl vs T* 



plunge, rn = 0, li = 2 ) 
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•Sure ( 4 . 12 ) C ’,1 vs T* ( plunge, in = 0 . h = 2 ) 
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Figure (4.13) 



Wake pattern ( plunge. 



m = 0.20833, 



h 
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Figure (4.14) 



Cl vs T* ( plunge, ui 



0.20833 , li = 2 ) 
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Figure (4.15) Cd vs T* ( m= 0.20833 , h = 
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Figure (4.16) Wake pattern ( plunge, in 



0.25 , li = 2 ) 
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Figure (4.17) 



Cl vs T* ( plunge, m = 0.25 , h 
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Non-Dimensional Time (T*) 




Figure (4.18) 



C<l vs T* 



( plunge, in = 0.25 



h =2) 
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Figure (4.19) Wake pattern (plunge, ni = 0.5 



h=2) 
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Figure (4.20) Cl vs T* ( plunge, in = 0.5 , li =2 
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Figure (4.21) 



Cd vs T* ( plunge, m = 0.5 



li = 2 ) 
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Figure (4.22) Average Propulsive Force Due to Plunge Versus Wake Phasing, in 
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Similar results are obtained for the pitch case. The airfoils studied are again 
NACA 0007 foils, pitching about the quarter-chord point with an amplitude of 
I degree and a reduced frequency of 0.0617. Figures (4 23) to (4.32) display 
the same phasing values as in the plunge case. These figures and the summa- 
rizing Figure (4.33) show again that values of ni near zero and one produce 
increased drag compared to the single airfoil case, but values of m near 0.5 
produce thaist, whereas the single airfoil produces drag. Both analyses, the anal- 
ysis by Wood and Couch [40] for the special case of two wakes only and the 
panel code, produce similar trends However, the quantitative agreement is not 
as good as in the plunge case. The precise reason is not sufficiently understood 
at this time It appears that the resolution of the suction peaks over the leading 
edge and therefore the computation of the thrust is more sensitive in the pitch 
case than in the plunge case. Furthermore, it must be remembered that Loewy's 
analysis for an infinite number of wakes and that by Wood and Couch [40] for 
two wakes are based on the flat-plate assumption, whereas the panel code requires 
a minimum blade thickness in order to obtain accurate results. A further evalu- 
ation of the panel code can be done for the drag produced by a single airfoil 
pitching about the leading edge. In this case the analytical theory of Garrick [5] 
and the computational results of Bosch [ 15 ] are available. It can be seen from 
Figure (4.34) that the trends are again in agreement, but the code results deviate 
from Garrick’s and Bosch’s results with increasing frequency. These deviations 
must be attributed to the geometry differences ( flat plate versus NACA 0007 ) 
and the higher amplitude used in the panel code. Note that the pitching airfoil 
starts to develop thrust only at a reduced frequency greater than about 1.8. 
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Figure ( 4.23 ) Wake pattern ( pitch, in = 0 , h -2 ) 
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Cl vs T* ( kl = k2 = 0.0617, a G = 1 0 , Wake Spacin 
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Figure (4.24) Cl vs T* ( pitch, m = 0 , li =2 ) 
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Non-Dimensional Time (T*) 
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Figure (4.25) Cd vs T* ( pitch, in = 0 , It =2 ) 
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Non-Dimensional Time (T*) 



Cl vs I * ( kl - k2 - 0.0617, a G -1°, Wake Spacing, h =2, m = 0.20833 ) 
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Figure (4.26) Cl vs T* ( 



pitch, in 



0.20833 , li =2 ) 
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Figure (4.27) (’d V s T* 



pitch, in = 0.2083 , Ii =2 ) 
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Cl vs T* ( kl = k2 = 0.0617, a Q = 1 0 , Wake Spacing, h =2, m = 0.25 ) 




* 

H 

v 

B 

H 

C3 

c 

o 

00 

a 

E 

O 

1 

c 

o 

z 



Figure (4.28) Cl vs T* 



pitch, in = 0.25 , li -2 ) 
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Figure (4.29) Cd vs T* ( pitch, m = 0.25 , h = 2 ) 
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Non-Dimensional Time (T*) 




Figure (4.30) Wake pnHern ( pitch, in - 0.5 , h —2 ) 
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Figure (4.31) (4 vs I* 



pilch, m = 0.5 . Ii -2 ) 
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Ctl vs T* ( kl = k2 = 0.0617, a Q = 1° , Wake Spacing, h 




Figure (4.32) Cd vs T* 



pitch, in = 0.5 , h -2 ) 
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I'igurc (4.33) Average drag vs wake phasing ( pitch only ) 
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Cduv Cdav vs k (Pitch, ea— O.AOA =5. 3. Sinclc Airfoil.) 
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Figure (4.34) Average drag versus reduced frequency ( single airfoil ) 
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Reduced Frequency (Theodorsen) 



V. PROPULSIVE EFFICIENCY 



AND FLUTTER ANALYSIS 

A. INTRODUCTION 

In the previous two chapters it was shown that a single airfoil which 
performs harmonic plunge oscillations generates a certain amount of forward 
thrust ( Katzmavr effect ) and that this thrust may be enhanced by the oscillatory 
wake shed from the preceding blade if the frequency ratio or phasing m is in 
the range 0.2 < m > 0.7 approximately 0.5. It was also shown that a similar 
enhancement occurs for pitching blades. However, pitching airfoils develop thrust 
only for relatively high values of reduced frequency, k > 0.6 by Garrick 
results [5], 

Furthermore, it is well known that a pitching airfoil may develop negative 
aerodynamic damping and thus experience single-degree-of-freedom flutter. As 
nown by Loewy [4], the possibility of flutter is enhanced by wake interference. 
This is an effect known as wake-induced flutter. It is therefore of interest to 
study the effect of various parameters on propulsive efficiency and on flutter 
instability in more detail. 
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B. PROPULSIVE EFFICIENCY 



The propulsive efficiency is defined as the ratio of the average work of 
propulsion to the work required to maintain the oscillation. If this factor is zero 
or negative, we expect energy to be extracted from the air and fed to the oscil- 
lating system. For a fiat plate performing a flapping motion given bv hit), where 
hit) = h (> e 1C01 the result for the horizontal force per unit span, averaged over 
one cycle, is given bv [42] : 

D = -TtpbcoV; F 2 (k) +G 2 (k) (5.1) 



where Flk) and Glk) are the real and imaginary parts of the Theodorsen lift 
deficiency function C(k). For the case of pitching motion, the average horizontal 
force per unit span length is given by the following equation, (5.2): 



— 3 -> 2 I 

D = -7tpb'co“a' 

0 



[r 



+ G“ 



1 

+ 

k^ 



- a 



+ 



-F 



F G 



where a G is pitch amplitude and (a) is the position of axis of rotation. a = - 1 
corresponds to the leading edge, a = + 1 is the trailing edge. The propulsive 
efficiency then is given by : 



n = 



DU 

W 



(5.3) 



W is the average work required to maintain the oscillation. For pitch, it is 
given by [42]: 



W = Ttpb^cCUcT { - 

o ' 9 








(5.4) 
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I in plunge, it is given bv 



W - apb'lwo ' * F ( k ) 
° k 



7tp b h“(D “UF (k) 



I ben t he efficiency factor. q, can be expressed as. for plunge 



n 



h 



l'^(k) i (i“(k) 
F (k) 



and for pitch 



(5.5) 



(5.6) 



■ 




Figure (5.1) shows the propulsive efficiency of a single plunging airfoil as 
a function of reduced fiequcncy Slowly flapping airfoils reach very high values 
of efficiency, approaching one. As the frequency of flapping is increased, the 
efficiency falls to 50 percent bitching airfoils, on the other hand, start out with 
an efficiency of no more than 50 percent at high frequency of oscillation which 
quickly falls to zero and to negative values with decreasing frequency. A negative 
propulsive efficiency means that no work is required to maintain the oscillation, 
hut instead energy is extracted from the airstream. This implies a loss of damping 
and thus the possibility of singlc-dcgrec-of-frecdom flutter in pitch. 
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Figure (5.1) Propulsive KITiciency n versus 1/k for a pltiiigini> airfoil 
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Figure (5.2) shows the propulsive efficiency of an airfoil which pitches about 
the leading edge The efficiency quickly drops to negative values as the fre- 
quency is reduced. As seen in Figure (5.2), infinitely large negative values are 
reached as the inverse frequency approaches 24, corresponding to a reduced fre- 
quency of 0.038. This is the flutter condition, as will be explained in more 
detail in the next section. 

Similar trends are found for other pitch axis locations, as seen in Figure 
5.3 ( pitch axis is at the quarter chord and mid-chord point ). The three quarter 
chord point is an exception. Here the propulsive efficiency is minus one at very 
high reduced frequency, which falls to larger negative values with decreasing 
frequency ( Figure 5.4 ). However, it is interesting to note that a small shift of 
the pitch axis of only 0.5 % in either direction restores the efficiency to 0.5 
for large frequencies, from which it starts to drop with decreasing frequency. 

These results show that there is a fundamental difference between airfoil 
plunge and pitch motions. Plunge oscillations of a single airfoil always produce 
thrust, and the efficiency decreases as the frequency is increased. Pitch oscil- 
lations, on the other hand, produce thrust only at relatively high values of re- 
duced frequency, k > 0 . 6 . This possibility of single-degree-of-freedom in pitch 
flutter is discussed in more detail in the next section. 
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I'iguie (5.2) I’ropulsivo Efficiency i] versus I /It for an airfoil pitching; 
nlioul (lie leading edge 
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Figure (5.3) Propulsive Efficiency i] versus 1/k for an airfoil pitching 
about the quarter - chord point and mid -chord point 
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Figure (5.4) Propulsive Efficiency r| versus l/k for an airfoil 
pitching about the three - quarter- chord point 
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C. SINGLE-DEGREE-OF-FREEDOM AIRFOIL FLUTTER 



Consider an airfoil hinged at its leading edge but elastically restrained from 
rotating about this axis by a torsion spring with constant K a (ft-lb /rad). The 
airfoil is placed in a low-speed airstream so that the unstrained position of the 
spring corresponds to zero angle of attack a. Figure (5.5). The equation of 
motion for this single-degree-of-freedom system is given by 



I a + K a = M . 

(t (X V 



(5.8) 



M v is the aerodynamic moment due to a (t) and I (X is the moment of inertia 
about the leading edge. Assuming that a(t)= a () e lon . then equation (5.8) can 
be written as 



a 





f CO ^ 


1 - 
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v co y 



0 ^ 



+ m = 0. 
v 



(5.9) 



where co a is the natural frequency of torsional vibration and is given by : 



0 ) = 




(5.10) 



and m v is shorthand for the dimensionless aerodynamic coefficient and can be 
written as [41]: 
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M 
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ttpb o> a e 
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icot a 
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(5.11) 
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where " a " is the elastic axis position and at the leading edge a = - I m v 
is a complex number and it is a function only of reduced frequency, k “ (ob/U. 



I'.quation (5.9 ) can be split into real and imaginary parts. 



Re { m } - Re 



0 ) ' 
<< 



Viipb 



0) J 



1m { m v | ~ 0 



(5.12 a) 
(5.12b) 



Flutter occurs at that value of the reduced frequency where the imaginary 
( out-of-phase) part of the aerodynamic moment becomes zero, provided that the 
corresponding real (in-phase) part yields a non-imaginary flutter frequency. 




Figure (5.5) An airfoil restrained to rotate about its leading edge 
in two-dimensional flow. 
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Figure 5.6 shows the variation of the real and imaginary ( damping ) pans of 
the pitching moment for an airfoil pitching about the leading edge. It is seen 
that the pitch damping goes through zero at a reduced frequency of 0 038. This 
value agrees with the previously found value using propulsive efficiency (energy 
extraction ) considerations. 

Figure 5.7 displays the comparison between the pitch damping coefficients 
(note that the pitch damping coefficient is defined by Theodorsen [3] as : 

CL 71 k f 1 A I 

c = M - — a i ( L +M. )+ -a: L, (5.13) 

m 2 « V2 / a v> J h . 

The difference in Cm shown in Figure 4.8 is due to a difference between the 
M a , L a , M h and terms in Equation 5.13 as determined by panel code 
( Riester [43] ) and Theodorsen’s theory [3], The imaginary part in Figure 5.7 
is obtained by Equation 5.11 computed by the panel code and Theodorsen’s 
analysis. It is seen that the two computations are in reasonable agreement, 
recognizing the difference in geometry ( NACA 0007 versus flat plate ). 

Systematic variation of the pitch axis location shows that zero pitch damping 
is possible for axis locations upstream of the quarter -chord point. No flutter is 
possible, on the other hand, for axis locations at or downstream of the quarter 
chord point. Figure 5.8. Loewy [4] has shown that wake interference greatly 
increases the possibility of flutter. The unsteady panel code provides an oppor- 
tunity to investigate and to extend this well-known finding. This is a problem 
which should be studied in more detail. 
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Real 8 Imag. pirt of m y vs k at x = - 1 




Reduced frequency k 



liiiinc (5.6) \ nr in I ion willi reduced frc(jiieney k of t lie real and imaginary 
parts of (lie non-dimensional aerodynamic moment m due to pitching of 
an aii foil about its leading edge in an incompressible flow. 
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C m vs k (pitching airfoil, a Q 




figure (5.7) ( oinpnrison of Cm (tiling.) vs Kp for NACA0007 Oscillating 
in pitch about t lie leading edge for a ( , = 1° and N =200 panels. 
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0 0.2 0.4 0.6 0.8 I 

Reduced frequency k 



I 1 it’ll rc (5.8) linag.(iny) vs k (for an airfoil pitching at ca = -0.5, 0. 0.0.5, and 1.0) 
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VI. RESULTS AND CONCLUSIONS 



A DISCUSSION OF RESULTS 

This research was undertaken, to study the effect of unsteady aerodynamics 
on reductions of main rotor shaft torque and engine power measured when 
the HHC system was applied to the OH-6A rotor. In Chapter III we study pro- 
pulsive force resulting from pure plunge oscillation of an airfoil in inviscid, 
incompressible tlow. The results of flat plate theory. Theodorsen [3] and Garrick 
[5], are used to validate panel code results. Comparison with Garrick showed 
good agreement in pure plunge, especially, for lower values of reduced frequency 
and amplitudes of oscillation. 

Figures 3.5 to 3 , show the average drag coefficient versus plunging 
amplitude for three values of reduced frequency, k = 0.1, 1. and 2 respectively. 

It is seen from Figure 3.5. ( reduced frequency, k is 0.1 ) that agreement 
between the panel code and Garrick [5] is very good for a range of non- 
dimensional plunging amplitude up to 15 to 20 % of the blade chord. On the 
other hand, for higher values of k, the two results are in good agreement for 
a range of plunging amplitude up to approximately 10%. Single airfoil code, was 
used to compute the time history of the drag (thrust) coefficient of the NACA 
0015 airfoil (the blade section of the OH-6A ). 
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Figure 3.12 and Table 3. 3 show the average drag coefficient versus reduced 
frequency for an NACA 0015 airfoil which executes plunging oscillation with 

1 " plunging amplitude. The OH-6A data presented in Table 3.1. Wood [42], 
and Figures 3 8 to Figure 3.11 from Ref [27], were used to estimate the 
propulsive force per blade obtained due to plunging oscillation for I P frequency 
case. In these figures, the difference between the deflections when the HHC 
is " on " and " off " represents the plunging amplitude. The controller phase is 
also significant Phase angle of 90° or 120° are those where helicopter vibra- 
tions were most sever, whereas at phase angle of 300° or 330°. the helicopter 
vibrations were significantly reduced. 

We look for the greatest power benefit to occur at the phase angle of 
lowest vibration, that is 300° or 330°. Table 3.2 shows the amount of propulsive 
force and torque generated at each blade segment. The results of this table 
showed that the contribution of IP (first harmonic) were 3 horsepower. Wood 
et al [27], applied the Garrick equation and considered the effect of blade plunge 
response in relation to the second to twelfth harmonic loads. He found for 2 P 
harmonic only (HHC "on", 330° phase at 60 knots ), the power gained was 
about 11.3 horse power. For more detail, the reader is encouraged to study Table 

2 of Ref. [27] 1 



1 Wood I:. R. . Higman. J. and Ramesh Kolar. " Higher Harmonic Control 
Promises Improved Dynamic Interlace Operations’ 1 AGARD Proceedings of 
the 78 th Flight Mechanics Panel, May 1991, Seville, Spain. 
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Chapter IV, considers the influence of the rotor wake interaction and examines 
whereas such interaction might cause favorable effects. First, the case of an iso- 
lated airfoil ( single airfoil oscillating in plunge) with absence of wake was 
considered ( or the wake located at great distance from rotor ). This is the "Katz- 
mayr effect" case where the resulting drag force is a propulsive force. It was 
the primary case considered by Garrick [5], 

The advantage of the Panel Code is that it provides us with time histories 
of unsteady lift, drag, and the wake trajectory by which we can better understand 
and analyze the problem. In each case, the wake pattern was plotted in order 
to be sure that the wake produced by the second airfoil would not interfere with 
the reference airfoil. Also the lift and drag (or the propulsive force) were plotted 
versus non-dimensional time, T* ( T* = t.U/c ). 

The closed form solutions are based upon a very thin flat plate airfoil 
undergoing very small motions when we model it in the unsteady panel code. 
For the calculations an NACA0007 airfoil was selected. The airfoil sets at zero 
mean angle of attack and the reduced frequency, k, is = 0.0617 

For the OH-6A, the value of reduced frequency , k at 70 % radial station , 
was calculated and found to be .1234 for the 32 Hz exciting frequency. The 
reference length used by the panel code is the blade chord or 2 b, so the 0.1234 
reduced frequency corresponds to 0.0617 per the panel code definition. Non- 
dimensional wake spacing, for single airfoil analysis, was h = 200 ( mathemat- 
ically equivalent to infinite spacing, or we say that the interaction between the 
two airfoil is very small) . 
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The amplitude ot' plunge oscillation was h 0 = 0.14 ( by the panel code 

■07, lor the same reason ). The panel code results are compared to other refer- 
ences in order to know the code limits. The comparison of the code results to 
that obtained by Ref. [44], is shown in Figure (4 2). This figure shows the 
relative lift coefficient of an airfoil ( a wing of relatively long span ) as it ap- 
proaches the ground for different values of angles of attack. The figure, showed 
that the results estimated by the panel code and the results obtained by Ref 
[44], are very close. For h/c * 0.5 the results obtained by the panel code 
are suspected. 

Another comparison is shown in Figure (4.4). This figure and Table 4 1 
show the correction factor, B, versus the chord-gap ratio in a biplane configura- 
tion. The correction factor, as defined by Ref. [45], is the reduction of the lift 
coefficient of an unstaggered biplane compared with that of monoplane at the 
same angle of incidence. On Figure (4.3) and Figure (4.5), the local speed ratio 
and pressure coefficient vs X/C are plotted and showed very good agreement 
with VVegley [46] and Bagley (1959), respectively. 

The panel code is suspect whenever the wake generated by an airfoil 
impinges on another airfoil. The numerical panel code method also encounters 
difficulty with discontinuities. For these reasons it is important to carefully check 
the time history plots after each computer run. 
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In Figure (4 8), the time history of the lift coefficient. Cj variation is 
plotted. Shown in Figure (4.9 is the drag coefficient, Cj versus the non- 
dimensional time T* We observe from Figure (4.8) that the lift varies as the 
frequency of oscillation, with a mean lift value of zero as we expect. Amplitude 
of oscillatory lift was found to be± 0.0510. Figure (4.9) clearly shows the " 
Katzmavr effect ". We noted that the calculated drag varies at twice the plunge 
oscillation frequency of the airfoil Also, observe that the airfoil is generating a 
propulsive force as indicated by the mean value of drag coefficient, which is 
seen to have a negative value. 

Looking to the computer printout at, at the non-dimensional time T* = 
0.0, steady state, the drag generated by the first airfoil was found to be 
Cds s =.000131 for NACA0007, which is the airfoil type that was used for 
most of our analyses at zero angle of attack , and a total number of panels 
N= 100 ( i. e 50 panels at the upper surface and 50 panels at the lower surface). 
For NACA0009, (this airfoil was used to compare the code results with that of 
Bosch’s for a pitching case), the steady state drag Cj ss was found to be .000196 
at N=100 and .000113 at number of panels N=150. 

It is known that, for an airfoil at zero angle of attack, in an inviscid 
incompressible flow field, steady state drag is zero, so in our calculations that 
steady state drag was subtracted from the drag generated at each time step. Shown 
in Figure (6.1) is the steady state drag coefficient vs angle of attack [37], Figure 
(6.1a) shows the variation of steady state drag coefficient vs angle of attack for 
three different airfoils; NACA0003, NACA0007, and NACA0012 at total number 
of panel of 200. 
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In this graph, it is shown that, for the thin airfoii NACA0003 at zero angle 
ot attack, the drag coefficient is very small but as the angle of attack increases, 
the drag coefficient highly increases more than in the cases of NACA0007 and 
NACA0012. Figure (6. lb) shows the effect of increasing the total number of pan- 
els N on the drag coefficient of the airfoil NACA0012. As the total number 
of panels increases, the steady state drag coefficient accumulated error decreases. 

Figures (4 10) through Figure (4 21), show the same case as above, that 
of an airfoil oscillating in plunge, but with the introduction of a wake layer 
in the near proximity of the reference airfoil (non-dimensionai wake spacing, h 
=2 ) for wake phasing, m = 0. 0.20833, 0.25. and m = 0.5. These correspond 
to phase angles of 0° , 75°, 90° and 180°, respectively. 

We see that the phase angles relate to the phase relationship between the 
reference airfoil and the wake layer of shed vorticitv. These phase angles were 
obtained in the panel code by establishing the proper time phase between the 
initiation of the wake vorticitv and the oscillation of the reference airfoii. This 
time phase is in turn a function of the frequency of oscillation which is directly 
related reduced frequency, k. For the results presented here, for k = .0617, this 
corresponds to 50.92 time units for m=0.0 (0° phase); 61.53 time units for m 
= 0.2083 ( 75° phase), and 76.38 time units for m = 0.50 ( 180° phase), etc. 
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Figure (6.1) Drag coefficient vs angle of attack 
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If the reduced frequency . k, changes, this corresponds to a change in 
frequency (at given chord length and free stream velocity ). For k =0.1234. this 
correspond 50.92 time units for m = 0.0, 56.23 time units for m = 0.2083, and 
63.65 time units for m = 0.5, etc. 

Figure (4 10) shows the wake position for m= 0.0, the case where the 
wake shed from the second airfoil at h = 2 is in phase with wake shed from 
the reference airfoil itself. In Figure! 4 21 ), the time history of lift variation 
is shown. 

This is indicated by the solid line for Cl j . The dashed line. Ch , indicates 
the variation of lift on a second airfoil located at the point of initiation of the 
wake layer at h = 2. Our discussion will focus airfoil], indicated by the solid 
line in the figures mentioned above. For an airfoil oscillating about zero angle 
of attack, the mean lift is zero. We found that, while the mean lift is zero, the 
amplitude of lift oscillation is smaller in the presence of the wake at zero 
phasing. From these figures, we conclude that the effect of wake vorticity at 
zero phasing is to reduce the " Katzmayr effect". 

Comparison of oscillatory lift values of the 75° wake phasing case, with 
the case of no wake or single airfoil, shows that they are identical, since the 
values at 75 degrees phase are the same as that for infinite wake. Table 6.1 
proves this result. 
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The time history of drag is seen to vary as twice the plunge oscillation 
frequency of the airfoil. We also observe that the airfoil is generating propulsive 
torce as indicated bv the mean value of drag, which is seen to have a negative 
value. Results of 75 degree wake phasing case, confirm that there is a phase 
angle, m, at which the effect of the wake vorticity on the reference airfoil gives 
the same results as when the wake is removed It has been shown at a phase 
angle of zero degrees that the shed wake vorticity increases the drag on the 
oscillating refeience airfoil, i.e. diminishes the "Katzmavr effect". 

Conversely we found that at a phase angle of about 180 degrees, there is a 
significant decrease in drag at the reference airfoil due to the wake vorticity. 
This implies an intermediate phase angle at which the drag is identical to the 
"Katzmavr effect". It turns out that such a phase angle exists. In fact, at this 
phase angle the lift, drag and pitching moment are found to be identical to the 
case of no wake, and match the "no wake" values identically. 

Table (6.1) show the aerodynamic coefficients Cl(l), Cl(2), Cm(l), Cm(2), 
Cd(l), and Cd(2) vs. Non-dimensional time. This table is a part of one of the 
output files. It shows these aerodynamic coefficient vs time from the nondimen- 
sional time unit 113.3 until 159.1. The phase angle at which this occurs is 75 
degree or m = 0.20833, for the case considered here. 
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fable (6.1) Aerodynamic coefficients vs Non-dimensional time (from time 
step 113 to 159) for the case of wake phasing, m =75 degree, wake spacing 
h = 2, reduced frequency k =0.0617 and plunging amplitude liliar = 0.14 ) 



TIME 



CM 1) 



CM2) CM ( 1 ) CM ( 2 ) CD ( 1 ) CD(2) 



111. 790Q.12 
! I 1 .563066 
l 1 5 .036000 
117.100731 
1 1 0 .302669 
1 19.655602 
1 20 .928535 
1 72.201169 
123.171103 
121.717337 
1 26 .07027 1 
127.293705 
128.566132 
129.030066 
131 . 112000 
132.381933 
133.657867 
1 31.930801 
1 36 .703735 
1 37.476669 
l 18.749603 
I 10.022537 
141 .295471 
142.560405 
113.041339 
145.114773 
146. 387207 
117.660141 
1 49.933075 
150 .206009 

151.170043 
1 52 .751877 
1 54 .02401 1 
155.297745 
166 . 570679 
1 57 .04 36 13 
159.116547 



0.006952 
-0.001073 
-0 .009081 
-0.016877 
-0.024266 
-0.031067 
-0.037111 
-0 .04274 7 
-0.016317 
-0.049307 
-0.051051 
-0.051534 
-0 . 050746 
-0.018701 
-0.045454 
-0.041083 
-0.035698 
-0.029433 
-0.022416 
-0.014910 
-0.007014 
0.001016 
0 . 009072 
0.016865 
0 . 024230 
0.030990 
0.036981 
0 .012057 
0.046093 
0.048996 
0.050693 
0.051145 
0.050343 
0 . 048306 
0.015081 
0.040757 
0.035427 



0.014965 
0.007130 
-0.000981 
-0.000871 
-0.016614 
-0.024005 
-0 .030777 
-0.036791 
-0.011898 
-0.045974 
-0.048918 
-0.050659 
-0.051152 
-0.050385 
-0.048379 
-0.045182 
-0.040872 
-0.035556 
-0.029364 
-0.022450 
-0.014902 
-0.007147 
0.000065 
0.008856 
0.016627 
0.023991 
0.030761 
0.036775 
0.041803 
0 . 045960 
0.048903 
0.050613 
0.051136 
0.050371 
0.040364 
0.045167 
0.040858 



-0.001387 
0.000665 
0.002701 
0.001671 
0 . 006534 
0 . 008236 
0.009738 
0.011001 
0.011994 
0.012693 
0.013077 
0.013140 
0.012878 
0.012296 
0.011412 
0.010245 
0.008825 
0.007188 
0.005375 
0.003431 
0.001404 
-0.000655 
-0.002696 
-0.004668 
-0.006521 
-0.008212 
-0.009700 
-0 .010948 
-0.011926 
-0.012610 
-0.012985 
-0.013039 
-0.012775 
-0.012196 
-0.011319 
-0.010163 
-0.008757 



-0.003432 
-0.001427 
0.000614 
0.002639 
0 . 004599 
0 . 006446 
0.008134 
0.009621 
0.010872 
0.011855 
0.012546 
0.012929 
0.012994 
0.012738 
0.012169 
0.011301 
0.010155 
0.008758 
0.007146 
0.005358 
0.003437 
0.001132 
-0.000608 
-0.002634 
-0.004593 
-0.006441 
-0.008128 
-0.009616 
-0.010867 
-0.011850 
-0.012541 
-0.012923 
-0.012988 
-0.012733 
-0.012164 
-0.011296 
-0.010150 



- 0.000011 
0.000000 
-0 . 000009 
-0 . 000036 
-0.000079 
-0.000133 
-0.000194 
-0.000255 
-0.000310 
-0.000355 
-0.000384 
-0.000391 
-0.000386 
-0.000358 
-0.000315 
-0.000260 
-0.000199 
-0.000138 
-0.000082 
-0.000038 
-0.000009 
0.000001 
-0.000007 
-0.000035 
-0.000078 
-0.000133 
-0.000194 
-0.000255 
-0.000310 
-0.000354 
-0.000382 
-0.000392 
-0.000384 
-0.000357 
-0.000314 
-0.000260 
-0.000199 



-0.000039 
- 0.000011 
0.000000 
-0.000008 
-0.000035 
-0.000078 
-0.000132 
-0.000193 
-0.000254 
-0.000309 
-0.000353 
-0.000382 
-0.000393 
-0.000384 
-0.000358 
-0.000315 
-0.000261 
- 0.000200 
-0.000139 
-0.000084 
-0.000040 
-0.00001 1 
0.000000 
-0.000008 
-0.000035 
-0.000078 
-0.000132 
-0.000193 
-0.000254 
-0.000309 
-0.000353 
-0.000382 
-0.000393 
-0.000384 
-0.000357 
-0.000315 
-0.000260 
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The proper phasing is achieved by positioning the reference airfoil 76.38 
time units to the rig. of the second airfoil, the airfoil that generates the wake. 
We found in the case of the plunging airfoil that the wake at zero degree phase 
increases the drag force, while the wake at 75 degree phase gives us close 
results to "no wake at all" case. So, we expect that the wake at 180 degree to 
increase the propulsive force. Shown in Figure (4.22), is the variation of the 
average propulsive force coefficient versus wake phasing, m. The circles in the 
graph represent panel code results for the parameters mentioned above. 

Superimposed on the plot of the figure are the analytical closed form 
results of the Loewy theory as modified for one wake, as shown in Figure 
(4.22) from Wood et al [38], Note the good agreement between the two sets 
of results. It is likely that the agreement can be further improved by increasing 
the number of panel from 100, used for these calculations, to 200 or 200 more. 

As was done previously for an airfoil oscillating in plunge, we will now 
consider the case of an airfoil oscillating in pitch with particular emphasis on 
drag or propulsive force. We will first look at the case without presence of 
layers of adjacent wake, then look at the wake effect including the important 
consideration of phasing. We will begin by considering the numerical results and 
comparing these results with Garrick, Bosch, for "no wake" case, and then com- 
pare with Loewy’ s theory when the wakes are included. 
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Considered is the case of an NACA0007 airfoil set at zero degree angle 
of attack The airfoil is oscillating with an amplitude 1 degree. The reduced 
frequency is taken at k = 0.061 7 and the non-dimensional wake spacing was 
h=200. which is numerically equivalent to having the wake at infinity'. The time 
history of the lift variation. Cl, and the time history of drag variation, Cd, reveal 
the following. We observe that the lift vanes as the frequency of oscillation, 
with a lift mean value of zero as we expect. 

In contrast to the lift time history, the drag time history indicates that the 
calculated drag varies at twice the frequency of the pitch oscillation The am- 
plitude is given by Cd= 0.000192. The mean drag value is found to be Cd - 
0.00009 and the mean lift is of course 0.0 at zero angle of attack. 

For an airfoil oscillating in pitch with the wake layer near the reference 
airfoil at non-dimensional wake spacing, h = 2. The case of zero degree phase, 
where the wake below the airfoil at h = 2 is in phase the wake shed from the 
reference airfoil itself. Review of the results indicates a reduction in the ampli- 
tude of lift variation for the m = 0.0 case when compared to the case with no 
wake at all (wake at infinity ). We see that the lift amplitude decreases. We also 
see an increase in mean steady drag accompanies the reduction in lift amplitude, 
for the case of zero wake spacing (m = 0). Proper phasing is achieved at m 
=0.5. At this value, we found that the lift enhancement in this case of pure 
pitch, is similar to what we found in the plunge case. 
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B CONCLUSIONS 



The most important conclusions that can be drawn from this research are; 

1 For drag, in both cases of pure plunge and pure pitch, wake phasing can 
increase or decrease the steady component of drag acting on the airfoil. The 

increase in the steady component of propulsive force due to wake phasing in the 

case of pure plunge was sufficient to significantly enhance the " Katzmayr " 
effect value of propulsive force. 

2. For an airfoil oscillating in plunge or pitch, for optimum reduction in drag 
and enhancement in vibratory lift, the phase angle of the wake vorticitv in the 
single wake case with respect to motion of the reference airfoil should be about 
180°. 

3. For an airfoil oscillating in plunge or pitch, the largest value of steady drag 

and smallest value of oscillatory lift occurs when the phase angle of the wake 

vorticity with respect to motion of the reference airfoil is about 0 °. 

4. For lift, in both plunge and pitch in the single wake case, wake phasing can 
increase or reduce the oscillatory lift acting on the airfoil. Where in the absence 
of wake vorticity, the effect of oscillations decreases the lift (lift deficiency ). 
With the wake present we also observe lift enhancement ( lift efficiency ). 

5. There exists a phase angle, 75° for the case, (reduced k = .0617, wake 
spacing, h = 2 ), where the effect of wake vorticity on the reference airfoil is 
identical to the case of the wake at infinity ( represented in the code by wake 
h = 200 ), or no wake vorticity at all. 
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Applying these findings to the OH-6A case, we can add the following 
conclusions 

6. Measured reductions in power are feasible when the "Katzmavr " effect and 
thesis - reported wake enhancements are included. 

7, Measured higher harmonic control results show least benefit at 100 knots 
due to the wake being transported greatest distance from rotor at this speed. 

8 Poor repeatability of measured open loop performance data is to be expected 
due to shifts in w'ake position 1 . 

9. For helicopters, it is standard rotor blade design to locate the center of 
gravity, aerodynamic center, and elastic axis at the quarter chord point to avoid 
flutter. This was illustrated in the analysis. Results of this thesis showed that 
flutter is impossible for this condition without wake interference. However, 
Loewy's results introduce the possibility of wake induced flutter. Future investi- 
gations, therefore, will have to explore this possibility in more detail. 



1 E. R. Wood. Max F. Platzcr. Ahmed Abourahma, Mark A. Couch 
On the Unsteady Aerodynamics of Higher Harmonic Control. 
Paper No. C 17, Nineteenth European Rotorcrat’t Forum, Ccmobbio 
(Comoi. Italy. September 1993. 
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C RECOMMENDATIONS 



1. A new flight test program, directed at studying helicopter performance, is 
needed. The NASA-Army OH-6A program was dedicated to vibration testing. 

2. Continued full scale HHC tests in the NASA Ames 40 ' x 80 ' wind tunnel is 
encouraged. In addition, it is suggested that tests be initiated in small university 
wind tunnels to verify the findings of this thesis. 

3. Controlled whirl tower testing is also recommended, since this allows careful 
measurement and control of performance parameters. 

4. A wind tunnel program is needed to explore similar cases and to verify the 
present research findings. 

5. The panel method is a significant analysis tool that can be applied to study and 
solve similar problems. It provides an excellent method for studying single or 
two-degree-of-freedom flutter problems. 
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